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{Aa}    elastic  incremental  stresses  in  the  "initial 
stress"  method. 

{Aa"}    initial  stresses  to  be  supported  by  equivalent 
nodal  forces 

{6}    nodal  displacements. 

e     uniaxial  strain. 

{c}    vector  of  strains  at  a  point. 

z  ultimate  uniaxial  compressive  strain  of  a  material, 

cu 

{e  }    vector  of  initial  strains. 

n  .    local  y-coordinate  of  node  ' i'. 

1  * 

y     coefficient  of  friction  between  brick  (block)  and 
mortar. 

v    Poisson's  ratio  of  a  material. 

£  ■    local  x-coordinate  of  node  'i1. 

i 

a     uniaxial  stress  at  a  point. 
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10}    vector  of  stresses  at  a  point. 
a,o     uniaxial  yield  stress. 

y 

a  normal  stress  at  interface  between  block  and 

n 

mortar. 

a     ultimate  strength  of  concrete  plate  in  uniaxial 
0  .       • 

compression,  psi. 

{0  }    vector  of  residual  stresses. 

a  peak  stress  in  biaxial  compression,  psi. 
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tensile  strength  of  a  material. 

shear  stress;  shear  stress  at  the  interface  between 
block  and  mortar. 

limiting  bond  shear  strength  of  mortar. 
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The  main  objective  of  this  investigation  is  to  obtain 
empirical  equations  to  predict  the  shear  capacity  of  nonre- 
inforced  concrete  masonry  elements.  An  analytical  investi- 
gation using  a  nonlinear  finite  element  analysis  and  an  ex- 
perimental investigation  were  carried  out  for  this  purpose. 

First,  an  experimental  program  was  initiated  to  obtain 
equations  for  predicting  the  strength  of  mortar  joints  under 
combined  compression,  shear,  and  bending.   The  test  speci- 
mens were  standard  three-block  prisms  with  three-eighths-inch 
mortar  joints.   Precompression  up  to  300  psi  was  applied. 
Based  on  test  results,  a  circular  interaction  curve  is  pro- 
posed for  predicting  the  strength  of  mortar  joints  under 
combined  stresses. 


In  the  second  phase  of  the  experimental  program,  the 
shear  strength  of  concrete  masonry  piers  made  with  one- 
fourth-scale  model  blocks  was  investigated.   The  main  vari- 
ables were  the  strengths  of  grout  and  mortar  and  the  magni- 
tude of  external  precompression.   A  relationship  between 
the  shearing  strength  of  a  grouted  pier  and  the  correspond- 
ing prism  strength  was  established.   Equations  are  proposed 
for  predicting  the  shearing  strength  of  a  nongrouted 
concrete  masonry  pier. 

A  nonlinear  finite  element  analysis,  using  the 
isoparametric  approach  and  higher  order  rectangular  elements, 
was  developed  to  predict  the  behavior  and  collapse  of  con- 
crete masonry  elements.   An  incremental  loading  procedure 
was  adopted  and  failure  was  investigated  at  integration 
points  due  to  yielding,  cracking,  or  crushing  of  the  masonry 
component.   Also,  criteria  were  developed  to  predict  the 
debonding  of  the  mortar  joints.   The  "initial  stress" 
approach  was  adopted  to  redistribute  the  released  stresses, 
due  to  cracking,  to  the  surrounding  elements.   The  initial 
stiffness  was  used  throughout  each  iterative  cycle  as  well 
as  through  all  the  load  increments.   Predictions  of  failure 
patterns,  collapse  load,  and  deformations  are  compared  to 
experimental  results. 


CHAPTER   1 
INTRODUCTION 

1.1   History :   Masonry  construction  is  perhaps  the  oldest 
building  system  employed  by  man  as  evidenced  by  the  historic 
remains  of  ancient  temples  in  India  and  Egypt.   These  early 
works  evolved  more  from  art  than  science  and  are  character- 
ized by  their  massiveness  and  quality  craftsmanship.   This 
traditional  masonry  construction  prevailed  throughout  the 
United  States  in  the  nineteenth  century  and  was  culminated 
by  the  construction  of  the  Monadnock  Building  in  Chicago, 
Illinois,  completed  in  1891.   This  16-story  skyscraper  has 
exterior  bearing  walls  varying  in  thickness  from  12  inches 
at  the  top  story  to  6  feet  at  the  base. 

Late  in  the  nineteenth  century,  the  concept  of  the 
structural  frame  in  structural  steel  or  reinforced  concrete 
replaced  masonry  bearing  wall  construction.   Masonry  became 
primarily  an  architectural  product  used  in  nonbearing  parti- 
tions and  as  a  facing. 

Expanded  research  in  building  methods  and  materials 
early  in  the  twentieth  century,  particularly  with  respect  to 
structural  concrete,  has  resulted  in  renewed  interest  in 
load-bearing  masonry.   Although  a  considerable  amount  of 
masonry  research  has  occurred  since  19  32  (4  3) ,  there 


exists  little  correlation  among  the  various  studies  conducted 
by  governmental,  promotional,  and  university  research 
agencies.   Each  study  has  of  economic  necessity  been  con- 
strained within  narrow  bounds  placed  on  the  variables.   For 
this  reason,  recommendations  based  upon  this  research  have 
been  purposely  conservative.   Much  reliance  has  been  placed 
upon  the  results  of  the  voluminous  research  and  experience 
with  concrete  and  reinforced  concrete;  and  research,  design, 
and  construction  procedures  for  reinforced  masonry  have  ad- 
vantageously followed  those  developed  for  reinforced  concrete. 
However,  there  are  important  differences  in  the  behavior  of 
masonry  and  concrete  that  must  be  taken  into  consideration. 

1.2   General  Remarks :   Concrete  masonry  walls  are  broadly 
classified  as  either  load-bearing  or  nonload-bearing  and 
further  described  in  various  ways  such  as  either  nonreinforced, 
solid  unit,  or  hollow  unit. 

The  general  behavior  of  high-rise, load-bearing  structures 
under  gravity  and  lateral  loads  is  the  combined  action  of 
floors,  bearing  walls,  and  shear  walls.   Floors  transmit  hori- 
zontal forces  by  diaphragm  action  from  the  exterior  walls  to 
the  shear  walls,  which  in  most  cases  are  also  the  bearing 
walls.   The  floor  system  must  be  sufficiently  rigid  to  serve 
as  a  diaphragm,  and  connections  must  be  adequate  to  transfer 
these  forces  to  the  shear  walls  which  carry  them  to  the 
foundation.   The  height  to  which  the  buildings  can  be 


constructed  depends  upon  the  strength  of  the  masonry 
materials,  the  spacing  of  intersecting  walls  and  floors,  their 
connection  to  each  other,  and  the  shape  of  the  structure. 

The  design  procedure  includes  an  investigation  of  the 
following  (41) (20) . 

1)  Bearing  capacity:   bearing  stresses  will  generally 
govern  required  block  strength  and  wall  thicknesses. 

2)  Stability  against  overturning:   overturning  resis- 
tance to  lateral  loads  depends  upon  the  shape  and  mass  of  the 
building.   Shear  walls  that  are  also  load-bearing  walls  are 
the  most  effective  structural  elements  for  developing  resis- 
tance to  overturning. 

3)  Shear  resistance  of  walls:   lateral  load  is  trans- 
mitted through  the  floors  to  those  shear  walls  parallel  to 
the  assumed  direction  of  the  lateral  load.   The  percentage  of 
total  lateral  load  carried  by  a  shear  wall  is  proportional  to 
this  "relative  rigidity"  with  respect  to  other  participating 
shear  walls. 

4)  Flexural  resistance  of  walls:   lateral  bending  of 
walls  can  be  produced  by  wind  loads  on  exterior  walls,  by 
eccentricity  of  loading,  and  by  insufficiently  rigid  floor 
diaphragms  and  shear  walls. 

5)  Floor-wall  connections:   finally,  because  the 
strength  and  stability  of  the  high-rise  building  depend  upon 
the  interaction  of  the  connecting  floor  and  wall  elements, 
connections  must  be  adequate. 


As  mentioned  previously,  shear  walls  are  designed  to 
resist  the  effects  of  lateral  forces  acting  on  buildings. 
The  lateral  forces  are  primarily  due  to  wind  or  earthquake. 
The  performance  requirements  for  shear  walls  under  wind 
loads  are  different  than  that  for  earthquakes  (3) .   Walls 
designed  for  wind  forces  have  to  meet  both  strength  and 
stiffness  requirements.   Walls  designed  for  earthquakes  must 
also  satisfy  requirements  of  ductility  and  energy  absorption, 
damping  characteristics  and  damage  control,  during  several 
cycles  of  inelastic  deformation  (1) . 

The  behavior  of  shear  walls  is  complicated  by  the 
influence  of  boundary  elements  and  multiple  openings  (14) (36) 
Figures  1.1  and  1.2  [taken  from  reference  (3)]  present  some 
typical  examples.   Lateral  loads  are  usually  introduced  into 
shear  walls  through  floor  slabs  framing  into  either  one  side 
or  both  sides  of  the  walls.   As  a  result,  the  lateral  loads 
tend  to  be  distributed  across  the  width  of  the  wall.   Trans- 
verse walls  or  columns  are  often  located  at  the  extreme 
edges  of  the  walls.   They  act  with  the  wall,  and  usually 
contain  most  of  the  flexural  reinforcement  resisting  the 
moment  due  to  the  lateral  forces. 

When  a  wall  contains  large  openings,  it  can  be 
considered  to  be  made  up  of  a  system  of  piers  and  spandrels. 
Each  individual  pier  or  spandrel  is,  in  effect,  a  shear  wall 
element,  with  a  shear  span  approximately  equal  to  one-half 
of  its  height  or  length,  respectively.   In  addition  to  shear, 


WALL  OR    COLUMN  — 


SHEAR  WAl! 


^%^%/^ 


Fig.  1.1.   Shear  wall  with  boundary 
elements . 
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Fig.  1.2.   Elements  of  shear  wall 
with  openings. 


piers  will  also  generally  have  tension  or  compression  caused 
by  gravity  and  overturning  forces  as  well  as  shrinkage, 
creep,  and  differential  settlement. 

It  is  noted  by  the  Joint  ASCE-ACI  Task  Committee  4  26 
on  Shear  and  Diagonal  Tension  that  the  shear  strength  of  a 
wall  is  of  interest  only  for  shear  span  to  depth  ratios  less 
than  2,  or  for  walls  with  a  flanged  cross-section  (3) . 

1. 3   Previous  Investigations: 

1.3.1   Strength  of  Mortar  Joints:   compared  with  the  vast 
number  of  tests  reported  on  concentrically  loaded  walls  with 
the  load  applied  vertically,  little  is  known  about  the 
strength  of  masonry  walls  with  the  load  applied  at  different 
inclinations  to  the  horizontal  joints.   However,  some  test 
results  are  reported.   Benjamin  and  Williams  (5)  carried  out 
tests  on  shear  couplets  of  two  bricks  bound  together  with  a 
mortar  joint.   Three  different  mortar  types  were  tested  with 
watered,  stiff  mud,  side-cut,  vacuum- treated  clay  bricks. 
The  test  results  showed  little  or  no  influence  of  brick  and 
mortar  compressive  strengths  on  the  couplet  bond  strengths 
in  tension  and  shear.   The  test  results  showed  that  the 
relationship  between  the  shear  stress,  fg  (shear  force 
divided  by  wall  area),  and  the  normal  stress,  f n ,  could  be 
expressed  in  the  form: 


fs  =  c  +  u  fn  (i.D 


where  \i    is  the  coefficient  of  friction  between  mortar  and 
brick  and  C  is  the  initial  bond  strength,  a  constant. 
Benjamin  and  Williams  obtained  a  value  of  C  =  150  psi  and 
y=  0.73.   They  also  concluded  that,  up  to  a  compressive 
stress  of  approximately  6  50  psi,  the  shear  strength 
increases  with  the  compressive  stress.   For  higher  compres- 
sive stresses,  the  apparent  shear  stress  still  increases, 
although  the  joint  has  already  failed  in  shear.   The  addi- 
tional strength  is  claimed  to  be  due  to  friction. 

Zelger  (44)  has  reported  tests  on  masonry  specimens  of 
the  type  shown  in  Figure  1.3.   Zelger  obtained  C  =  2  kg/sq  cm 
and  \i   =  0.5  in  his  tests.   Yorulmaz  and  Sozen  (44)  tested 
masonry  specimens  of  model  bricks  0.53"  x  0.86"  x  1.87"  in 
size.   The  results  obtained  from  test  specimens  of  the  type 
in  Figure  1.3  gave  C  =  150  psi  and  y  =  0.46. 

Haller  (22)  adopting  the  test  setup  shown  in  Figure  1.4, 
arrived  at  the  following  empirical  relationships: 


f   =  35/f   +  280  -  540  psi  (1.2) 

s       n  c 


f   <  200  psi 

n  —     c 


for  normal  quality  masonry  consisting  of  cored  bricks  (3300 
psi  to  6500  psi)  and  cement- lime-sand  mortar  (1225  psi) . 

If  the  Haller  formula  is  approximated  to  linear 
relationship,  the  following  equation  is  obtained: 
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Fig.  1.3.   Type  of  test  specimens  used  by 
Zelger. 
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Fig.  1.4.   Test  specimen  adopted  by  Haller. 


f   =  50  +  0.88  f  (1.3) 

s  n 


f   <  200  psi 

n  — 


The  tests  of  the  type  shown  in  Figure  1.3  give  the 
most  representative  values  for  shear  strengths  of  the  mortar 
joints  in  masonry,  since  disturbances  caused  by  the  testing 
machine  platens,  etc.,  are  much  less  likely  to  occur  in 
this  type  of  test  compared  with  the  couplet  type  tests. 

From  the  limited  number  of  tests  mentioned,  it  seems 
reasonable  to  assume  that  the  bond  or  shear  failure  of  a 
mortar  joint  in  brick  masonry  follows  Equation  (1.1)  in  a 
range  of  approximately  2%  to  15%  of  the  compressive  strength 
of  the  masonry.   C  is  of  the  order  of  2%   to  3%  of  the  com- 
pressive strength  of  masonry.   However,  since  the  couplet 
tests  indicate  that  the  compressive  strength  of  the  mortar 
has  no  influence  on  the  bond  strength,  the  shear  strength 
should  be  tested  when  high  stresses  are  employed. 

For  compressive  stresses  lower  than  approximately  2% 
of  the  compressive  strength  of  the  masonry,  and  for  pure 
tensile  stresses,  the  shear  strength  falls  below  that  calcu- 
lated from  Equation  (1.1),  as  revealed  by  couplet  tests  and 
model  masonry  tests.   The  pure  tensile  bond  strength  is 
greatly  influenced  by  workmanship  and  wetness  of  the  brick?. 
A  suction  rate  of  20  g/min  or  less  seems  to  give  maximum 
bond,  although  saturated  bricks  produce  close  to  maximum  bond. 
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For  high  compressive  stresses,  the  apparent  shear  strength 
again  is  lower  than  calculated  from  Equation  (1.1). 

Hedstrom  (24)  reports  load  tests  of  concrete  masonry 
walls  with  constant  wall  dimensions  but  with  mortar  bed 
joints  in  90°,  45° ,  and  0°  inclination  to  the  axial  load, 
which  was  applied  parallel  to  the  longer  side  of  the  walls. 
The  tensile  bond  strength  obtained  with  the  two  types  of 
mortar  was  tested  on  masonry  prisms  of  two  blocks  in 
bending.   A  plot  of  bond  shear  strength  vs.  compressive 
strength  yielded  C  =  48.5  psi  and  y  =  0.84  for  M-type 
mortar  and  C  =  24  psi  and  y  =  0.92  for  S-type  mortar.   The 
figures  are  supported  by  too  few  tests  to  be  conclusive. 

1.3.2   Racking  Tests:   the  present  standard  racking  test 
described  in  ASTM  E  72-68,  Method  of  Conducting  Strength 
Tests  of  Panels  for  Building  Construction,  provides  only  a 
relative  measure  of  the  shearing  or  diagonal  tensile 
resistance  of  a  wall.   Results  of  this  test  are  consequently 
valid  only  for  comparison  purposes  and  are  not  suggested  for 
determination  of  design  values. 

In  this  method  of  test  (Figure  1.5),  horizontal 
movement  of  the  wall  specimen  (81  x  8'),  due  to  the  horizon- 
tal racking  load  at  the  top  of  one  end,  is  prevented  by  a 
stop  block  at  the  bottom  of  the  other  end.   To  counteract 
rotation  of  the  specimen  due  to  this  overturning  couple, 
tie  rods  are  used  near  the  loaded  edge  of  the  wall  specimen. 
Under  racking  load  these  rods  superimpose  an  indeterminate 
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Fig.  1.5.   Arrangement  of  racking  test  and  the  force 
distribution  on  the  specimen. 
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Fig.  1.6.   Test  for  diagonal  tensile  strength  of 
brickwork  and  stress  distribution. 
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compressive  force  which  suppresses  the  critical  diagonal 
tensile  stresses  and  increases  the  load  required  to  rack 
the  specimen.   A  typical  mode  of  failure  is  indicated. 
The  obtained  apparent  shear  strength  as  calculated  from 
f  =  T/ld  is  usually  of  the  order  of  25  psi  to  300  psi 
under  laboratory  conditions. 

For  concrete  masonry  walls,  the  racking  strength  was 
reported  by  Fishburn  (18)  to  be  25  psi  to  50  psi  for 
masonry  walls  having  a  compressive  strength  of  390  psi  to 
470  psi  giving  a  racking  strength  of  about  7%  to  10%  of 
compressive  strength. 

1.3.3   Circular  Shear  Specimens  (28):   in  these  tests 
(Figure  1.6),  a  15"  diameter  specimen  is  tested  in  compres- 
sion with  the  line  of  load  at  45°  to  the  bed  joints.   As 
shown  in  Figure  1.6,  the  diametrical  stresses  are  largely 
tensile  over  the  central  80%  of  the  specimen.   The  tensile 
stress  is  approximately  constant  for  about  60%  of  the  dia- 
meter and  may  be  calculated  by  the  following  equation: 


f.  =  |£r  (1.4) 

t     TTDt 


where  P  equals  load  at  rupture,  in  pounds,  D  equals  the 
diameter  of  the  specimen,  in  inches,  and  t  equals  the  thick- 
ness of  the  specimen,  in  inches. 
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1.3.4   Square  Shear  Specimens  (6)  (16):   in  examining  the 
damage  done  during  earthquakes,  it  was  noted  that  cracks  in 
shear  walls  were  frequently  of  a  diagonal  nature,  so  fre- 
quently that  they  were  called  typical  "X  cracks"  (Figure  1.7) 
With  this  as  a  starting  point,  the  theory  was  early  recog- 
nized that  the  force  of  quake  working  against  the  static 
resistance  of  a  building  would  produce  a  racking  effect; 
this  in  turn  would  be  resisted  by  the  diagonal  strength  of 
the  wall  member.   It  was  found  that  the  "X"  cracking  devel- 
oped in  a  diagonal  tension  test. 

Photoelastic  analysis  through  the  use  of  polarized 
stress  panels  was  used  to  demonstrate  the  validity  of  this 
theory  (Figure  1.8).   The  stress  distribution  in  the  wall 
panel  was  shown  to  be  the  same  as  in  the  diagonal  test 
panel.   It  was  then  decided  that  the  proper  test  would 
consist  of  breaking,  by  diagonal  loading,  4'  x  41  panels 
incorporating  desired  variables  (Figure  1.9).   In  this 
method,  the  test  results  are  susceptible  to  stress  analysis. 
In  addition,  they  are  more  reproducible  and  thus  more 
reliable  for  comparison  and  design  data  purposes. 

The  square  specimen  is  placed  in  the  testing  frame  so 
as  to  be  loaded  in  compression  along  a  diagonal,  thus  pro- 
ducing a  diagonal  tension  failure  with  the  specimen  split- 
ting apart  along  the  loaded  diagonal.   The  shear  strength, 


14 


.Jl j _Ji 

_  jc  zr 


Fig.  1.7.   Typical  X-cracks  in  a  wall  damaged 
by  an  earthquake. 


Fig.  1.8.   Fringe  patterns  obtained  in  a  photoelastic 
analysis  of  a  model  of  a  wall  with  an 
opening  and  a  square  panel  under  diagonal 
compression. 
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Fig.  1.9.   Diagonal  tension  test  on  a 
square  panel. 
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Fig.  1.10.   Test  setup  (schematic)  adopted 
by  Meli  and  Reyes  for  testing 
mortar  joints  under 
compression  and  shear. 
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v  ',  is  determined  from  the  equation  (6) (27) 


O^TOp  (1#5) 


where  F  equals  the  diagonal  compressive  force  or  load,  in 
pounds,  t  equals  the  thickness  of  wall  specimen,  in  inches, 
and  1  equals  the  length  of  a  side  of  a  square  specimen,  in 
inches. 

1.3.5   Horizontally  and  Vertically  Loaded  Wall  Without  Frame: 
the  load  carrying  capacity  of  a  wall  subjected  to  a  horizon- 
tal load  at  one  of  the  upper  corners  is  governed  mainly  by 
the  shear  and  tensile  strength  of  the  bedjoints  at  the 
foundation  of  the  wall  (44) .   By  precompression,  for  example, 
by  dead  load  from  slabs  and  walls  above,  the  strength  is  in- 
creased in  a  manner  similar  to  that  described  for  masonry 
specimens  loaded  with  an  inclined  load. 

Murthy  and  Hendry  (44)  report  "1/6-modal"  tests  on 
three  bay,  one-story,  shear  walls,  0.669"  thick,  about  16" 
in  height  and  length.   The  bricks  had  an  average  strength 
of  4421  psi  and  the  cement- lime- sand  mortar  about  1200  psi. 
The  horizontal  shear  strength  was  tested  for  various  addi- 
tional loads  up  to  180  psi,  and  the  following  relationship 
was  established: 


f   =  30  +  0.5  f  (1-6) 

s  n 


f   <  180  psi 

n  —     c 
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Benjamin  and  Williams  (5)  tested  model  walls,  without 
frames  and  without  vertical  loads,  and  found  apparent  shear 
strengths  of  15  psi  to  30  psi. 

1.3.6   Test  on  Small  Masonry  Assemblages:   at  the  University 
of  Mexico,  Meli  and  Reyes  (39)  conducted  tests  on  small 
assemblages  for  investigating  the  mechanical  properties  of 
masonry.   Three  tests  were  found  to  be  the  most  useful:   a 
small  prism  subjected  to  axial  compression,  a  wallette  under 
diagonal  compression,  and  a  three  unit  assemblage  subjected 
to  shearing  of  the  joints.   Results  of  the  prism  test  were 
related  with  the  behavior  of  masonry  walls  under  vertical 
loads.   The  remaining  two  tests  were  related  with  the 
behavior  of  walls  under  lateral  loads.   Tests  were  performed 
on  a  large  number  of  specimens  built  with  commonly  used 
types  of  masonry  units  and  mortars. 

Based  on  test  results,  it  was  found  that  the  prism  with 
height/thickness  ratio  of  4.0  gave  satisfactory  and  uniform 
index  to  the  resistance  to  axial  load  of  masonry. 

Figure  1.10  shows  the  schematic  test  setup  adopted  for 
the  shearing  tests  on  joints  with  precompression.   The 
results  were  expressed  in  the  form  of  Equation  (1.1)  with 
C  =  1.8  kg/sq  cm  and  y  =  0.8  for  concrete  blocks  with  dif- 
ferent types  of  mortar  whose  strengths  varied  from  151  kg/ 
sq  cm  to  43  kg/sq  cm.   Coefficient  of  friction  was  found  to 
be  a  very  uniform  property  for  the  different  types  of  bricks 
and  concrete  blocks  adopted.   The  values  were  approximately 
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0.7  in  all  the  cases.   However,  the  value  of  adhesion  varied 
depending  on  the  type  of  brick  and  mortar.   At  low  levels  of 
confinement,  the  results  varied;  however,  uniformity  in  re- 
sults was  obtained  with  high  level  of  precompression. 

Figure  1.11  shows  the  different  types  of  specimens 
adopted  for  the  diagonal  compression  tests.   For  each  series 
of  specimens  shown  in  Figure  1.11,  four  different  types  of 
mortars  and  seven  types  of  masonry  units  were  used.   For 
each  combination  of  materials,  there  were  three  specimens 
of  each  type.   The  object  of  the  series  of  the  tests  was  to 
study  the  effect  of  the  variation  of  the  height/width  ratio 
(h/1)  of  the  panel  and  number  of  joints  in  the  specimen  on 
the  resistance  of  the  assemblage.   The  results  showed  that 
each  type  of  specimen  had  a  definite  mode  of  failure  inde- 
pendent of  the  type  of  mortar  used.   In  general,  the  failure 
occurred  by  shearing  along  the  joints  for  long  specimens  and 
by  diagonal  tension  for  specimens  with  high  h/1  ratio.   The 
type  of  failure  was  not  always  perfectly  defined.   In  many 
cases,  the  crack  crossed  the  joints  and  the  blocks  partially. 
Qualitatively,  it  could  be  said  that  when  the  failure  was  by 
diagonal  tension,  the  resistance  was  relatively  uniform  for 
similar  specimens  and  did  not  depend  much  on  the  type  of 
mortar  used.   On  the  contrary,  when  the  failure  was  by  shear, 
the  dispersion  of  results  was  very  high.   It  was  found  that 
the  resistance  increased  with  increase  in  h/1  ratio. 
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However,  the  effect  of  this  increment  was  of  very  little 
influence  when  the  failure  was  by  diagonal  tension  and 
very  much  noticeable  when  the  failure  was  by  shearing  along 
the  joints.   It  might  be  due  to  the  fact  that  the  h/1  ratio 
controls  the  value  of  compression  normal  to  the  joints  that 
causes  the  effect  of  friction  to  be  developed. 

All  the  specimens  considered  in  this  investigation 
were  nongrouted.   It  was  found  that  the  strength  in  diagonal 
tension  was  equal  to  the  square  root  of  the  compressive 
strength  of  the  prism  when  the  stresses  were  expressed  in 
kg/sq  cm. 

For  the  shear  failure  along  the  joints,  a  relation 
similar  to  that  of  Polyakov  (42)  was  proposed: 

f   =  0.8/(1  -  0.9  u  h/1)  d-7) 

s 

1.3.7   Reinforced  Concrete  Masonry  Walls  in  Shear:   Figure 
1.12  shows  the  steel  arrangement  of  reinforced  concrete 
masonry  walls  tested  by  Converse  (12) .   The  basic  mode  of 
failure  was  one  of  diagonal  tension.   Walls  under  group  'B' 
showed  an  increase  in  strength  of  38%  over  group  'A'.   It 
suggested  that  the  increase  in  strength  was  nearly  propor- 
tional to  the  areas  of  additional  steel,  irrespective  of 
position.    Due  to  the  difference  in  number  of  bars,  no 
direct  comparison  could  be  made  of  the  relative  effect  of 
vertical  and  horizontal  steel,  but  indications  were  that 
they  were  equally  effective. 
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Scrivener  (48)  confirmed  the  previous  findings  in  his 
tests.   The  objectives  of  his  tests  were  to  determine: 
1)  the  pattern  of  behavior  as  the  percentage  of  reinforcing 
was  increased,  2)  the  relative  effectiveness  of  vertical 
and  horizontal  reinforcing,  and  3)  the  difference  between 
the  behavior  of  walls  where  the  vertical  steel  was  periph- 
eral and  where  the  steel  was  distributed  over  the  length  of 
the  wall.   The  following  conclusions  were  drawn  from  these 
series  of  tests: 

1)  Vertical  and  horizontal  reinforcing  are  equally 
effective  in  producing  satisfactory  crack  behavior  and 
failure  loads. 

2)  Walls  with  evenly  distributed  reinforcing  have  a 
later  onset  of  severe  cracking  than  walls  where  the  rein- 
forcement is  concentrated  in  the  wall  periphery. 

3)  With  a  low  percentage  of  reinforcing,  failure 
occurs  soon  after  the  onset  of  severe  cracking.   With  higher 
percentages  of  reinforcing,  the  failure  load  is  much  greater 
than  the  load  causing  severe  cracking. 

4)  Higher  failure  loads  were  obtained  with  walls  with 
higher  percentages  of  reinforcing  up  to  0.3%  of  the  gross 
cross-sectional  area  of  the  wall.   Above  this  percentage, 
additional  reinforcing  had  little  effect  on  the  failure 
load.   From  the  walls  with  the  optimum  (0.3%)  or  higher  per- 
centage of  reinforcing,  the  ultimate  horizontal  shear  stress 
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(ultimate  load  divided  by  the  gross  cross-sectional  area 
of  the  wall)  was  found  to  be  170  psi. 

Confirmation  of  the  last  conclusion  can  be  had  from 
the  test  results  of  Schneider  (4  5) ,  who  found  a  maximum 
effective  quantity  of  reinforcement  of  0.2%.   In  his  tests, 
a  racking  load  only  was  applied,  but  sufficient  peripheral 
vertical  reinforcement  was  placed  in  the  walls  to  prevent 
this  steel  being  stressed  beyond  its  yield  point  when  sub- 
jected to  the  tensile  forces  induced  by  the  maximum  over- 
turning moment.   Schneider's  walls  failed  in  shear  with 
the  typical  diagonal  cracking.   The  difference  between 
Schneider's  and  Scrivener's  test  walls  lies  in  the  boundary 
conditions  on  the  vertical  sides.   Schneider  also  found 
that,  workmanship  and  reinforcement  remaining  the  same,  the 
shear  resistance  of  stack  and  running  bond  and  stack  bond 
masonry  block  walls  was  about  the  same. 

1.3.8   Effect  of  Wall  Openings:   openings  in  shear  walls 
are  mainly  due  to  doors,  corridors  and  mechanical  duct 
space.   When  the  opening  is  relatively  small  and  spaced,  at 
least  a  distance  equal  to  the  size  of  the  opening  in  each 
direction,  its  influence  on  the  behavior  of  the  structure 
is  negligible  (29).   Figure  1.13  shows  an  elevation  of  a 
typical  shear  wall  in  an  apartment  building  using  an  8" 
concrete,  flat  plate  slab  construction.   The  opening  shown 
at  the  center  indicates  the  corridor  at  each  floor.   In 
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Fig.  1.13,   Continuous  opening  in  a 
shear  wall. 
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Fig.  1.14.   Staggered  opening  in  a 
shear  wall. 
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apartment  buildings,  the  opening  would  normally  be  from 
top  of  the  floor  to  bottom  of  the  floor  above.   In  apart- 
ment buildings  where  the  architectural  planning  frequently 
permits  the  shear  wall  to  extend  from  one  face  of  the 
building  to  the  other,  openings  such  as  shown  in  Figure 
1.13  limit  the  full  utilization  of  the  entire  shear  wall  as 
a  unit.   The  connecting  slab  at  each  opening  is  relatively 
very  flexible.   As  a  result,  the  shear  wall  acts  as  two 
individual  shear  walls  on  either  side  of  the  corridor 
opening.   Although  for  medium  height  buildings,  this  inef- 
ficiency does  not  seriously  affect  the  economy  of  the 
entire  structure,  for  heights  above  40  stories,  its  effect 
on  overall  economy  becomes  significant.   As  a  solution  to 
this  problem,  Khan  (29)  proposes  staggering  of  such  open- 
ings at  alternate  floors  in  order  to  maintain  the  struc- 
tural continuity  of  the  entire  shear  wall.   Figure  1.14 
illustrates  the  proposed  arrangement. 

In  office  buildings  the  requirement  for  mechanical 
duct  space  under  the  floor  slab  makes  the  problem  of  open- 
ings different  from  that  in  the  apartment  buildings.   A 
hung  ceiling  is  almost  always  necessary.   The  door  opening 
size  generally  allows  a  2  ft  to  5  ft  connecting  beam  over 
each  opening.   The  proper  analysis  and  design  of  the  con- 
necting beam  is  important  because  the  beam  not  only  con- 
nects the  adjacent  parts  of  a  shear  wall  for  monolithic 
action,  but  also  redistributes  loads  in  different  parts  of 
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the  shear  wall.   Giri javallabhan  (21),  on  the  basis  of  an 
analysis  using  the  Finite  Element  Method,  found  that  one 
of  the  most  influencing  factors  on  the  distribution  of 
stresses  and  deformations  in  the  shear  wall  was  the  stiff- 
ness of  the  lintel  beam.   He  varied  the  depth  of  the  lin- 
tel beam  and  studied  the  influence  of  the  variation  on  the 
overall  behavior  of  the  shear  wall. 

Kokinopoulos  (30)  conducted  a  photoelastic  analysis  on 
models  of  single  story  walls  with  openings  to  study  the 
effect  of  size  of  opening  on  the  stress  distribution  in 
walls.   Schneider  (46)  conducted  full  scale  tests  on  piers 
in  a  shear  wall  with  openings.   This  is  the  only  experi- 
mental investigation  on  this  problem  which  utilized  full 
scale  tests.   This  investigation  is  described  in  the  next 
subsection. 

1.3.9   Shear  in  Concrete  Masonry  Piers:   this  investigation 
was  carried  out  to  estimate  the  capacity  of  concrete  masonry 
piers,  functioning  within  the  confines  of  a  shear  wall  in  a 
building,  to  resist  lateral  load  effects.   To  simulate  as 
nearly  as  possible  conditions  that  occur  in  an  actual  wall, 
such  as  degree  of  restraint,  amount  of  reinforcing,  manner 
of  vertical  load  imposition,  and  magnitude  of  secondary 
stresses,  the  pier  was  considered  along  with  the  wall  around 
it  (Figure  1.15) .   The  fully  restrained  configuration  was 
either  10'  8"  or  11'  4"  high.   The  cantilever  pier  was  either 
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Fig.  1.15.  Forces  acting  on  a  pier 
element  in  a  wall  with 
opening. 
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DIAGONAL  LOAD  FRAME 
(PIER  DI-TAILS  VARY) 


Fig.  1.16.  Test  setup  adopted  by 

Schneider  for  full  scale 
tests  on  piers. 
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7'  4"  high  and  8'  0"  wide  or  10'  0"  high  and  10'  8"  wide. 
Also  tested  was  a  set  of  4  8"  square  concrete  masonry  assem- 
blies.  A  diagonal  loading  frame  was  used  to  apply  the 
loading  (Figure  1.16).   Pipe  columns  were  used  to  maintain 
the  geometry  of  the  openings. 

The  important  variables  considered  were: 

1)  a/D  ratio,  where  a  is  the  distance  from  the  point 
of  inflection  to  the  spandrel  restraint,  and  D,  the  overall 
pier  depth. 

2)  Amount  of  web  reinforcement  (both  horizontal  and 
vertical) . 

3)  Amount  of  jamb  reinforcement. 

4)  External  axial  compressive  stress. 

5)  Nongrouted  panel  behavior. 

The  average  strength  of  blocks  used  varied  from  13  38 
psi  to  2962  psi.   The  28-day  mortar  strength  varied  from 
2487  psi  to  5116  psi.   The  strength  of  grout  prisms  varied 
from  1789  psi  to  3414  psi.   Reinforcing  steel  conforming 
to  ASTM  A  15,  had  a  yield  strength  of  55,000  psi  and  ulti- 
mate tensile  strength  of  80,000  psi. 

The  following  are  the  main  findings  of  this  investigation 

1)  The  shear  strength  increased  with  a  decrease  in  the 
a/D  ratio,  and  this  rate  of  increase  jumped  sharply  below  an 
a/D  ratio  of  0.5:1. 

The  very  consistency  of  the  test  results  throughout  the 
range  of  a/D  ratios  selected  for  analysis  suggested  the 
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following  relationships  for  average  ultimate  shear  stress 
of  a  pier  contained  within  a  shear  wall,  where  web  rein- 
forcement is  not  provided: 

ULTIMATE  MASONRY  SHEAR  STRESS, 
v  ,  psi 
M/VD  (a/D) 2 

Fixed  Pier  Elements: 

0.10:1  <  M/VD  <  0.50:1  V/tD  =  310  -  350  M/VD 

0.50:1  <  M/VD  <  1.50:1  V/tD  =  152.5  -  35  M/VD 

1.50:1  _<  M/VD  V/tD  =  100  psi 

Cantilever  Pier  Elements: 

1.0:1  <    M/VD  <  3:1  V/tD  =95-15  M/VD 

3.0:1  <  M/VD  V/tD  =  50  psi 

where  M  equals  spandrel  moment,  V  equals   base  shear,  D 
equals  overall  pier  depth,  a  equals  distance  from  the  point 
of  inflection  in  the  span  of  the  pier  to  the  fixed  end,  and 
t  equals  net  thickness  of  wall. 

2)   The  presence  of  adequate  horizontal  web 
reinforcement  materially  increased  the  shear  resistance  of 

the  pier. 

The  following  relationships  between  the  M/VD  ratio  and 
ultimate  shear  stress  were  proposed  for  piers  with  horizontal 
web  reinforcement. 
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ULTIMATE  MASONRY  SHEAR  STRESS 
WITH  ADEQUATE  WEB  REINFORCEMENT, 

M/VD  (a/D)  Vuf  PS1 

Fixed  Pier  Elements: 

0.10:1  <_   M/VD  <  0.50:1  V/tD  =  347.5  -  225  M/VD 

0.50:1  _<  M/VD  <  1.50:1  V/tD  =  290  -  110  M/VD 

1.50:1  <  M/VD  <  2.00:1  V/tD  =  200  -  50  M/VD 

2.00:1  <  M/VD  V/tD  =  100  psi 


3)  Vertical  steel  did  not  seem  to  function  effectively 
as  web  reinforcement. 

4)  Assuming  that  enough  jamb  steel  was  present  to 
resist  the  end  moments,  any  further  increase  did  not  alter 
the  pier  resistance  appreciably. 

5)  The  existence  of  a  bedjoint  fracture  at  a  foundation 
did  not  seem  to  impair  the  ability  of  the  panel  to  resist 
lateral  loads.   However,  if  a  bedjoint  crack  occurred  at  the 
center  of  a  square  panel,  where  tensile  stresses  are  maximum, 
its  shear  resistance  was  drastically  reduced. 

6)  The  energy  absorbing  ability  of  an  adequately  rein- 
forced masonry  pier  was  well  demonstrated.   As  the  shape  of 
the  load  deflection  curves  indicated,  these  piers  were  able 
to  absorb  a  great  deal  of  inelastic  strain  energy  without 
collapsing  or  even  spalling  seriously. 

7)  Concrete  masonry,  if  properly  reinforced,  exhibits 
a  tendency  toward  a  ductile  behavior  throughout  the  loading 
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sequence.   It  can  sustain  significant  proportions  of  the 
ultimate  load  well  into  the  inelastic  regions  (beyond  the 
first  significant  crack)  while  undergoing  rather  large 
lateral  deflections.   It  also  exhibits  effective  dampening 
characteristics,  especially  after  cracking. 

8)   On  the  basis  of  defining  ductility  as  the  ratio 
of  the  total  deflection  experienced  to  the  deflection  at 
the  first  shear  crack,  which  was  assumed  to  be  the  incep- 
tion of  inelastic  deformation,  the  ductility  factor 
exceeded  two,  which  is  considered  a  desirable  level. 
1.3.10   Strength  of  Masonry  under  Combined  Compression  and 
Shear:   Sinha  and  Hendry  (50)  propose  that  brickwork  sub- 
jected to  combined  compression  and  shear  exhibits  two  dis- 
tinct types  of  failure: 

1)  Shear  failure  at  the  brick/mortar  interface.   The 
shear  strength  consists  of  initial  bond  shear  and  the  resis- 
tance, proportional  to  the  normal  stress,  due  to  friction 
between  brick  and  mortar. 

2)  Diagonal  tensile  cracking  through  bricks  and  mortar 
governed  by  constant  maximum  tensile  stress  or  strain. 

On  the  basis  of  tests  on  circular  shear  specimens, 
Sinha  and  Hendry  found  the  diagonal  tensile  strength  of 
brickwork  to  be: 


£fc  -  2.0  *- 
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Let  f  equal  precompression  applied  to  brickwork  and 

f   equal  shear  stress.   If  it  is  assumed  that  failure  is 
s 

determined  at  a  certain  stage  by  the  criterion  of  maximum 
tensile  stress,  then: 


f 
2  -  -^  =  constant  (1.9) 


ft  =  ^f  2/4  +  f 

T-       y  '  < 


For  failure: 


fs  1  »   fy 


It  is  assumed  that  the  above  condition  will  be  fulfilled  by 

two  values  of  f  :   f  i  and  fy2.   Between  the  precompressive 

stresses  f  i  and  f  2  failure  of  the  structure  will  occur  by 

attaining  maximum  tensile  strength.   Below  and  above  this 

range,  failure  will  be  governed  by  shear  at  the  brick/mortar 

interface.   Precompression  above  f  2  will  suppress  the 

inherent  failure  due  to  diagonal  tension  and  modify  its 

value.   At  very  high  precompression  values,  the  failure  of 

the  brickwork  will  take  place  in  compression. 

Since  f  =  C  +  V  f  ,,  where  C  is  the  initial  bond  shear 
s         y 

strength  between  brick  and  mortar, 

ft  =  /f  ,2/4  +  (C  +  u  f  i>'  -  fy*/*  (1"10) 
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When  f     =   jj   f   2 , 


ft   =    /fv iV4    +    (y    f    2)2    "    fya/2  (1.11) 


Knowing  f  ,  C  and  y,  f  i  and  f  2  can  be  calculated.   Thus, 

t         y      y 

the  ultimate  shear  strength  may  be  calculated  from  the 
following  formulae: 


fs  =  C  +  y  f  (1.12) 


for  y  f   <  f  1 ,  and, 
y  —  y 


ft  -  ^f  2/4  +  f  2  -  V2  (1'13) 


for  f  1  <  f   <  f  ,/  and, 

y1  —   y  —   y 


f  =  y  f  (1.14) 


for  f  2  £  f  £  compressive  strength  of  brickwork. 


However,  in  a  recent  study,  Smith  et  al.  (51)  have 
concluded  that  the  diagonal  tensile  strength  of  brickwork 
is  approximately  equal  to  the  tensile  strength  of  brickwork 
or  mortar,  whichever  is  weaker.   This  conclusion  was  based 
on  an  analysis  of  a  masonry  circular  shear  specimen  using 
the  Finite  Element  Method  supported  by  experimental  data. 
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1.4   Objectives  and  Scope  of  Present  Investigation;   A 
survey  of  the  available  literature  revealed  a  need  for  more 
elaborate  research  on  shear  strength  of  concrete  masonry 
walls  with  openings.   To  meet  this  objective,  an  experimen- 
tal and  an  analytical  investigation  is  attempted  in  this 
dissertation.   In  the  experimental  investigation,  it  was 
decided  to  adopt  the  test  specimen  proposed  by  Schneider  (46) . 
Thus,  the  main  objective  of  the  experimental  investigation 
was  to  obtain  empirical  equations  to  predict  the  shear  capa- 
city of  grouted  and  nongrouted  piers  without  reinforcement. 
Since  the  effect  of  various  types  of  configurations  and 
reinforcement  had  already  been  established,  it  was  decided 
to  restrict  the  investigation  to  one  particular  configuration. 
The  main  variables  considered  for  grouted  piers  were  the 
grout  and  mortar  strengths;  for  the  nongrouted  piers,  the 
effect  of  external  precompression  normal  to  bedjoints  was 
treated  as  the  main  test  parameter  for  different  types  of 
mortar  adopted. 

Since  the  strength  of  mortar  joints  is  a  main  factor 
determining  the  behavior  of  a  nongrouted  masonry  element, 
an  extensive  study  was  initiated  to  obtain  equations  for 
predicting  the  strength  of  mortar  joints  under  combined 
compression,  bending  and  shear.   This  investigation  is 
described  in  Chapter  3 . 
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A  full-scale  testing  of  piers  was  not  possible 
because  of  limitations  of  capacity  of  testing  machines, 
space,  and  cost.   The  present  investigation  on  piers  was 
restricted  to  1/4-size  models.   Chapters  4  and  5  describe 
the  model  tests  on  piers. 

While  investigating  a  complex  phenomenon  such  as  the 
behavior  of  masonry,  if  a  suitable  analytical  model  could 
be  proposed,  it  would  facilitate  understanding  the  stress 
distribution  in  the  structure  more  thoroughly.   Such  a 
model  was  attempted  using  the  Finite  Element  Method.   The 
model  chosen  was  capable  of  predicting  the  different  failure 
modes  associated  with  masonry  and  a  nonlinear  finite  element 
analysis  was  adopted  accordingly  to  determine  the  ultimate 
load  and  failure  pattern.   The  analytical  results  were  com- 
pared with  experimental  ones.   The  analytical  investigation 
is  described  in  Chapters  6  through  8. 


CHAPTER  2 
EQUIPMENT,  MATERIALS  AND  TESTING  TECHNIQUES 

2.1   Testing  Machines  and  Other  Equipment:   The  Civil 
Engineering  Laboratory  is  equipped  with  a  hydraulic  press 
of  300,000  lb  capacity  and  a  mechanical  press  of  160,000 
lb  capacity.   The  clearance  of  these  two  machines  allowed 
good  observations  of  all  sides  of  the  specimen  being 
tested.   For  testing  the  model  reinforcement,  the  10,000 
lb  capacity,  Instron  Machine,  Model  TTC ,  in  the  department 
of  Metallurgical  Engineering  was  used.   The  machine  was 
provided  with  suitable  gripping  devices  for  clamping  small 
test  specimens  and  an  automatic  recorder  for  plotting  the 
load-extension  curves. 

The  investigation  on  the  strength  of  mortar  joints 
was  carried  out  in  the  300,000  lb  hydraulic  press.   The 
higher  clearance  and  lower  load  ranges  required  for  the 
model  tests  determined  the  use  of  the  160,000  lb  mechanical 
press  for  testing  model  piers. 

Sieves,  mechanical  shakers,  and  balances  were  available 
for  analysis  of  sand. 

An  electric  rotary  mixer  of  5  cu  ft  capacity  was  used 
in  mixing  mortars  for  building  3-block  prisms  used  for  deter- 
mining the  strength  of  mortar  joints  under  combined  stresses. 
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The  same  mixer  was  also  used  for  mixing  cement- sand  mortar 
for  making  spandrels  for  piers.   A  standard  flow  table, 
standard  molds  for  mortar  cubes,  paper  molds  for  mortar 
cylinders,  and  different  size  rods  were  used  for  mortar  and 
grout  control. 

A  large  moist  room  was  used  for  storage  of  mortar 
cubes  and  cylinders,  and  sufficient  storage  room  was  avail- 
able for  the  masonry  prisms  and  model  specimens. 

A  6"  diameter  brass  disc  (normally  used  for  capping 
concrete  cylinders) ,  together  with  four  1/4"  x  1/4"  brass 
bars,  formed  a  mold  for  the  sulfur  caps  of  the  model  con- 
crete blocks  and  prisms. 

Dial  gages  with  a  least  count  of  0.001"  were  used  for 
measuring  deflection.   They  were  mounted  on  magnetic  stands 
for  easy  removal. 

The  electrical  strain  gages  used  were  BLH's  SR-4  fixed 
on  the  test  surface  with  Duco  cement.  The  strain  indicator 
was  a  portable  Baldwin  Type  N.  The  strain  gages  were  wired 
to  the  indicator  through  a  switch  selector. 

2.2   Concrete  Blocks:   The  concrete  blocks  used  for  the  first 
investigation  had  nominal  dimensions  of  8"  x  8"  x  16".   The 
net  area  of  a  block  based  on  the  average  of  five  specimens 
was  61.4  sq  in.   Selected  blocks  were  capped  with  sulfur  and 
tested  according  to  ASTM  C-90  and  C-140.   The  average  com- 
pressive strength  based  on  net  area  was  6500  psi.   The  ini- 
tial rate  of  absorption  (IRA)  of  the  block,  according  to 
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ASTM  C-67,  is  measured  as  the  amount  of  water  initially 
absorbed  by  a  dry  unit  when  it  is  partially  immersed  in 
water  to  a  depth  of  1/8"  for  a  period  of  one  minute.   IRA 
is  measured  in  grams  per  minute  per  30  square  inches.   The 
test  was  conducted  in  the  following  manner. 

Four  steel  bars,  1/4"  x  1/4"  in  cross-section  and  4" 
long,  were  provided  with  needles  1/8"  high.   A  steel  pan, 
which  area  was  much  larger  than  the  gross  area  of  the  con- 
crete block,  was  chosen  and  placed  on  a  level  surface.   The 
steel  bars  were  positioned  on  the  pan  in  such  a  manner  that 
the  concrete  block  could  rest  on  them.   Water  was  allowed 
to  stand  in  the  pan  till  the  needle  points  were  just 
immersed.   The  previously  weighed  block  was  placed  in  posi- 
tion over  the  steel  bars,  and  the  water  supply  was  continued 
to  cope  with  the  absorption  of  the  block.   After  one  minute, 
the  block  was  removed,  the  immersed  surface  wiped,  and  the 
block  reweighed.   The  initial  rate  of  absorption  of  the 
blocks  varied  between  12  g  to  17  g/min/30  sq  in. 

An  indirect  tension  test  was  conducted  to  find  the 
tensile  strength  of  the  block  (49)  .   Two  test  methods  were 
devised  for  splitting  hollow  concrete  block.   The  first 
method,  shown  in  Figures  2.1(a)  and  2.1(b),  was  arranged  so 
that  the  block  could  be  split  twice.   The  load  was  applied 
through  hexagonal  bars  with  3/4"  flats,  first  through  one 
cell  and  then  through  the  other  cell  of  the  block.   The 
second  method,  shown  in  Figure  2.1(c),  was  arranged  to  load 
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both  cells  simultaneously.   Round  bars  of  1-3/8"  diameter 
were  used  to  distribute  the  load  in  this  test.   Usually- 
only  one  split  occurred  as  is  shown  in  Figure  2.1(d).   The 
results  of  the  tests  are  given  in  reference  (49) .   The 
indirect  tensile  strength  of  the  block  was  calculated  using 
the  relation: 


2P 


where  P  is  the  splitting  load  and  A  the  sectional  area  of 
split.   In  the  range  of  the  block  strengths  tested,  it 
appears  that  the  split  tensile  strength  is  approximately 
five  times  the  square  root  of  the  compressive  strength  as 
determined  by  the  standard  block  compression  test.   For 
the  blocks  used  in  this  investigation,  the  indirect  tensile 
strength  was  found  to  be  405  psi. 

2.3  Mortar :   The  cements  and  sand  were  provided  by  local 
suppliers  in  Gainesville,  Florida.   The  cements  were  manu- 
factured by  Florida  Portland  Cement.   Portland  Cement  Type 
I  and  Masonry  Cement  were  used  in  the  mortar  mixes  in 
accordance  with  ASTM  C-27  0-6  3.   The  granulometry  of  the 
sand  is  shown  in  Figure  2.2.   For  the  first  investigation, 
two  types  of  mortars  were  used.   The  following  are  their 
proportions  by  volume. 
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Fig.  2.2.   Sieve  analysis  of  sand  for  masonry  mortar. 


41 


Type  I 
Type  II 


PC 

1 

1 


MC 

1 
1 


Sand 

4.5 

6 


An  initial  rate  of  flow  of  about  100%  was  adopted.   Two- 
inch  cubes  were  molded  and  tested  according  to  ASTM  C-270. 
The  average  compressive  strength  of  Type  I  mortar  was  1948 
psi  and  that  of  Type  II  917  psi. 

2.4  Model  Concrete  Blocks :   The  model  concrete  blocks  were 
provided  by  the  National  Concrete  Masonry  Association, 
Arlington,  Virginia.   They  were  modelled  to  be  one-fourth 
the  size  of  the  full  block  of  nominal  dimensions  8"  x  8" 
x  16".   Thus,  the  model  block  had  nominal  dimensions  of  2" 
x  2"  x  4".   A  typical  block  is  shown  in  Figure  2.3.   The 
average  length  of  the  block  was  3.9",  width  1.9",  and 
height  1.85".   The  net  area,  based  on  an  average  of  six 
specimens,  was  4.16  sq  in.   The  model  blocks,  capped  with 
0.15"  thick  sulfur  capping,  were  tested  in  axial  compression. 
The  average  of  thirteen  tests  yielded  a  compressive  strength 
of  2688  psi  based  on  net  area.   The  coefficient  of  variation 
was  20.15%. 

The  absorption  of  the  model  block  was  11.9  8%.   It  was 
difficult  to  measure  the  initial  rate  of  absorption  for  the 
model  blocks  because  1)  their  rate  of  absorption  was  quite 
high,  and  2)  a  suitable  depth  of  immersion  of  surface  could 
not  be  defined  for  models. 
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Fig.  2.3.   Model  concrete  block 


Fig.  2.4.   Model  reinforcement  (0.147"  dia.) 
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The  tensile  strength  of  the  block  could  not  be 
determined  by  an  indirect  tension  test  because  of  the  un- 
certainties involved  in  choosing  the  size  of  splitting 
bars  for  models.   Also,  a  suitable  testing  machine  was  not 
available.   A  value  equal  to  five  times  the  square  root  of 
the  compressive  strength  of  the  block  was  assumed. 

In  order  to  obtain  the  stress-strain  curve  in 
compression,  two  SR-4  gages  were  mounted  centrally,  one  on 
each  longitudinal  side.   The  gage  length  was  0.2".   The 
stress-strain  curve  obtained  is  shown  in  Figure  2.5. 

2.5  Model  Reinforcement:   Number  2  bars  of  40  ksi  grade 
steel  were  used  for  main  reinforcement  in  spandrel  beams  of 
model  piers.   For  vertical  and  shear  reinforcements,  0.147" 
diameter,  high  strength  Duro-Wall  bars  were  used.   A  typical 
bar  is  shown  in  Figure  2.4.   The  bar  had  minute  depressions 
0.3"  long  and  alternate  projections  0.1"  long.   The  yield 
and  ultimate  strengths  of  the  bar  were,  respectively,  68  ksi 
and  75  ksi.   A  typical  stress-strain  curve  is  shown  in 
Figure  2.6. 

2.6  Model  Mortar:   Based  on  trial  mixes,  the  proper  grade 
of  sand  for  modelling  Type  I  and  Type  II  mortars  was  arrived 
at.   The  sand,  of  gradation  shown  in  Figure  2.2,  was  sieved 
through  a  set  of  sieves  (nos.  8,  16,  30,  50,  100,  and  200, 
in  that  order) .   Trial  mixes  were  made  with  the  sands  re- 
tained on  sieves  nos.  50,  100,  and  200.   It  was  found 
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Fig.  2.5.   Stress-strain  curve  for  model  concrete  block. 
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Fig.  2.6.   Stress-strain  curve  for  model  reinforcement. 
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that  for  the  same  water-cement  ratio  and  proportions  of 
ingredients,  the  strengths  of  mortar  obtained,  using  the 
sand  retained  in  no.  100  sieve,  compared  favorably  with 
those  of  Type  I  or  Type  II,  as  the  case  may  be.   Hence, 
this  grade  was  chosen  as  the  model  sand.   No  attempt  was 
made  to  model  cements.   An  initial  flow  of  120%  was  adopted. 
Compressive  strength  of  the  mortar  was  determined  by  testing 
standard  2"  cubes.   The  tensile  strength  was  obtained  by 
splitting  cylinders  3"  in  diameter  and  6"  long.   The  model 
mortar  properties  are  summarized  in  Table  2.1.   Plates, 
5.75"  x  5.75"  x  0.625",  made  of  the  different  types  of 
mortar  and  cured  in  air,  were  tested  in  uniaxial  compression. 
The  plate  compression  strength  was  used  in  the  analytical 
investigation  described  in  Chapter  8.   The  stress-strain 
curves  for  the  two  types  of  mortars  obtained  on  the  basis  of 
tests  on  3"  x  6"  cylinders  are  given  in  Figure  2.7. 

2.7   Model  Grout  Mixes :   No  gravel  could  be  used  in  model 
grout  because  of  difficulties  in  obtaining  proper  compaction 
due  to  the  small  size  of  the  cells.   Hence,  the  model  grout 
consisted  of  only  cement  and  sand.   Standard  2"  cubes  molded 
in  standard  molds  and  cured  in  air  were  used  for  measuring 
the  compressive  strength  of  the  grout.   Due  to  the  uneven 
surfaces  obtained  for  cubes  molded  by  using  the  concrete 
blocks  as  molds,  the  test  results  on  those  cubes  were  not 
satisfactory.   Plates,  5.75"  x  5.75"  x  0.625",  were  tested 
in  uniaxial  compression.   Their  compressive  strengths  were 
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Fig.  2.7.   Stress-strain  curve  for  model  mortar. 
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TABLE  2.1 
PROPERTIES  OF  MODEL  MORTAR  AND  GROUT 


Initial 


28-day  Strength 


Type 

Proportions  by 
Volume 

Rate  of 
Flow 

% 

Plate 
psi 

Cube 

Comp. 

psi 

Split 
Cylr. 

PC 

MC 

Sand 

psi 

Mortar : 

MS 

1 

1 

4.5* 

120 

656 

1695 

200 

MW 

1 

1 

6   * 

120 

485 

905 

135 

Grout: 

GW 

1 

5   + 

450 

570 

75 

GM 

1 

4   + 

760 

800 

205 

GS 

1 

3   ** 

1000 

1850 

322 

Note:   *  Gradation  of  sand:  passing  through  sieve  no.  50  but 
retained  on  sieve  no.  100. 
+  Gradation  of  sand:  passing  through  sieve  no.  30  but 
retained  on  sieve  no.  50. 
**  Gradation  of  sand:  as  shown  in  Figure  2.2. 
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Fig.  2.8.   Stress-strain  curve  for  model  grout. 
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Fig.  2.9.   Stress-strain  curve  for  spandrel  mortar. 
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required  for  use  in  analysis  later.   The  indirect  tensile 
strength  was  obtained  by  splitting  cylinders  3"  in  diameter 
and  6"  long.   The  properties  of  grout  mixes  adopted  are 
summarized  in  Table  2.1.   For  grout  types  GW  and  GM,  the 
sand  used  was  of  such  a  grade  it  passed  through  sieve  no. 
30,  but  was  retained  on  sieve  no.  50.   For  type  GS ,  the 
sand  used  for  mortars  Type  I  and  Type  II  was  adopted. 
Initially,  the  required  amount  of  water  was  determined  by 
trial  for  each  type  so  that  the  grout  could  easily  be  poured 
in  the  cells  without  separation  of  its  components.   The 
stress-strain  curves  obtained  using  3"  x  6"  cylinders  are 
shown  in  Figure  2.8. 

2.8   Mortar  Mix  for  Spandrel:   A  mortar  mix  of  the  following 
proportions  by  weight  was  used  in  forming  the  spandrel  of 
the  pier  specimen  used  in  model  tests  described  in  Chapter  4: 
Portland  Cement: Sand  (river  sand  used  in  concrete  mixes)  =  1:3, 
A  water-cement  ratio  of  0.55  by  weight  was  adopted.   The 
compressive  strength  of  the  mix  was  3830  psi  and  the  tensile 
strength  400  psi.   Compression  tests  on  a  plate  of  5.75"  x 
5.75"  x  0.625"  yielded  2300  psi.   The  stress-strain  curve 
obtained  by  testing  3"  x  6"  cylinders  under  uniaxial  com- 
pression is  shown  in  Figure  2.9. 


CHAPTER  3 
STRENGTH  OF  MORTAR  JOINTS  UNDER  COMBINED  STRESSES 

3.1   Scope:   This  study  was  designed  to  investigate  the 
influence  of  precompression  upon  the  shear  and  flexural 
tensile  strengths  of  mortar  joints  in  concrete  block 
walls.   First,  the  strength  of  the  joint  is  studied  under 
a  state  of  pure  shear.   Secondly,  this  is  followed  by  a 
study  of  pure  bending.   In  the  third  phase,  the  strength 
of  the  joint  is  studied  under  combined  bending  and  shear 
for  a  known  amount  of  precompression.   Combining  the 
above,  the  study  is  directed  toward  obtaining  interaction 
diagrams  that  can  be  used  to  predict  behavior  under  com- 
bined loading. 

The  structural  bond  between  mortar  and  units  is  an 
important  factor  in  the  structural  behavior  of  masonry 
walls  where  the  masonry  is  subjected  to  forces  which  pro- 
duce tensile  and/or  shearing  stresses  in  the  joints. 
Factors  which  seem  to  influence  the  bond  of  masonry  mor- 
tars are  (13)  (25):  1)  composition,  2)  suction  rate  (IRA), 
3)  initial  moisture  content  of  the  block,  4)  compressive 
strength,  5)  air  content,  6)  initial  flow  of  mortar,  and 
7)  curing  of  specimens.   Shear  bond  may  be  considerably 
increased  by  the  presence  of  compressive  stress  normal  to 
the  shearing  force. 
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In  the  present  investigation,  only  the  mortar 
composition  and  strength  were  considered  as  variables. 
Other  factors  such  as  suction  rate,  initial  flow,  block 
characteristics  and  curing  procedures  were  standardized 
and  held  constant. 

3.2  Test  Specimens:   Test  specimens  were  standard  3-block 
prisms  with  3/8"  mortar  joints.   All  specimens  were  air- 
cured  under  laboratory  conditions.   The  initial  rate  of 
absorption  of  the  blocks  varied  between  12g  to  17  g/min/ 
30  sq  in.   The  net  area  of  a  block,  based  on  the  average 
of  five  specimens,  was  61.4  sq  in.   Since  the  blocks  were 
fully  bedded,  this  value  was  taken  as  the  area  of  a  mor- 
tar joint.   Type  M  mortar  was  used  with  two  compositions 
The  properties  of  mortars  used  are  summarized  in  Table  3.1. 

3.3  Strength  of  Mortar  Joints  under  Compression  and  Shear: 
A  simple  method  was  developed  for  applying  shear  to  the 
joints  while  a  constant  uniform  compressive  stress  was 
maintained  across  the  joints.   The  specimens  were  supported 
and  loaded  in  such  a  way  that  the  joints  were  at  the  points 
of  contraflexure.   The  testing  arrangement  adopted  by  Base 

(4)  was  used  with  certain  modifications.   The  test  setup 
adopted  is  shown  schematically  in  Figure  3.1  and  positioned 
in  the  hydraulic  press  in  Figure  3.2.   The  axial  force 
across  the  section  was  applied  through  steel  rods  and  was 
varied  by  tightening  nuts  at  the  threaded  ends.   The  amount 
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TABLE  3.1 
PROPERTIES  OF  MORTAR 


Initial    28-day  Strength 


Proportion  by 

Rate  of 

Cube 

Split 

Mortar 
Type 

Volume 

Flow 

% 

Comp. 
psi 

Cylr. 

Specimens 

PC 

MC 

Sand 

psi 

PSM  1-12 

II 

1 

1 

6 

90 

864 

172 

PSM  13-25 

II 

1 

1 

6 

103 

932 

198 

PSM  26-37 

I 

1 

1 

4.5 

96 

2138 

290 

PSM  39-41 

II 

1 

1 

6 

110 

956 

- 

PSM  42-44 

I 

1 

1 

4.5 

100 

1758 

— 

55 


Test    Spe 


1"- thick  bearing     plate 
3/4"    threaded   rods 


Celotex  Board 


MJJJ.tf/f  TrhfrrTT 


SHEAR    FORCE    DIAGRAM 


EENM  KG   MOMENT    DIAGRAM 


Fig.  3.1.   Setup  for  studying  the  strength  of 
mortar  under  pure  shear. 
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Fig.  3.2.   Test  setup  for  finding  the  strength  of  mortar  joints 
under  combined  compression  and  shear. 


Fig.  3.3.   Shear  failure  through  mortar  joint, 
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of  prestress  developed  in  the  rods  is  found  by  measuring 
the  strains  with  electrical  strain  gages.   On  each  rod, 
two  strain  gages  were  mounted  diametrically  opposite  to 
each  other  and  connected  in  series  to  compensate  for 
bending  effects. 

Table  3.2  presents  the  summary  of  test  results  for 
the  low  strength  mortar  and  Table  3.3  for  the  high  strength 
mortar.   For  the  low  strength  mortar,  at  precompressions 
up  to  220  psi,  bond  failures  were  observed;  but  at  a  pre- 
compression  of  300  psi,  shear  failure  through  the  mortar 
frequently  occurred.   Figure  3.3  shows  an  example  of  shear 
failure  through  the  joint.   However,  for  the  high  strength 
mortar  (within  the  range  of  precompression  adopted) ,  the 
failure  pattern  was  always  one  of  bond.   Testing  beyond 
the  range  of  precompression  (300  psi  for  the  high  strength 
mortar)  was  not  possible  due  to  flexural  failure  of  the 
concrete  block. 

Figure  3.4  shows  the  plot  between  precompression  and 
the  shearing  strength  of  mortar  joints.   It  is  evident 
that  the  effect  of  precompression  is  to  increase  the 
shearing  strength  of  the  joint  considerably.   It  can  also 
be  observed  that  a  higher  strength  mortar  does  not  appre- 
ciably increase  the  shear  strength.   This  statement  cannot 
be  generalized  because,  within  the  range  of  precompression 
applied  thus  far  in  this  investigation,  bond  failures  were 
usually  observed.   It  is  quite  possible  to  have  an  increased 
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Fig.  3.4.   Effect  of  precompresslon  on  the  shearing 
strength  of  mortar  joints. 
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shear  strength  due  to  increase  in  strength  of  mortar  at 
higher  precompressions ,  because  the  failure  would  likely 
be  due  to  shear  failure  of  the  mortar. 

Using  the  method  of  least  squares,  the  best-fit  curve 
for  the  observed  data  was  found  to  plot  as  straight  lines 
with  the  following  equations : 

fso  =  °'606  fa  +  34  (3.1) 

for  1:1:6  mortar,  and, 

f    =  0.641  f   +  54  (3  2) 

so  a  \j . &) 

for  1:1:4.5  mortar. 

The  above  equations  may  be  used  for  predicting  the 
shearing  strength  of  mortar  joints  with  a  precompression 
up  to  300  psi  for  Type  M  mortar  when  blocks  having  an  IRA 
of  12  g  to  17  g/min/30  sq  in  are  used. 

The  specifications  (2)  limit  the  allowable  shearing 
stress  to  34  psi  for  Type  M  mortar,  irrespective  of  the 
strength  of  mortar  or  level  of  precompression.   Test  re- 
sults to  date  indicate  clearly  that  the  specifications  are 
conservative. 
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3. 4   Strength  of  Mortar  Joints  under  Compression  and 
Bending;   The  test  setup  shown  in  Figure  3.1  was  modified 
to  produce  bending  in  the  specimen  with  flexural  stresses 
acting  normal  to  the  bedjoints.   The  load  corresponding 
to  the  first  crack  was  taken  as  the  failure  load  since 
beyond  the  initial  failure,  the  load  could  have  been  shared 
by  the  prestressing  rods.   Table  3.4  presents  the  results 
of  the  tests.   The  maximum  bending  moment  corresponding 
to  the  failure  load,  when  divided  by  the  section  modulus 
of  the  joint,  yielded  the  value  of  modulus  of  rupture 
shown  in  Table  3.4.   Adopting  the  cross-section  of  mortar 
joint  shown  in  Figure  3.5,  the  value  of  section  modulus 
corresponding  to  the  extreme  fiber  was  obtained  as  117  cu  in. 
The  results  are  plotted  in  Figure  3.6.   Failure  at  the  joint 
was  precipitated  by  a  break  in  bond  between  the  mortar  and 
the  blocks. 

The  modulus  of  rupture  appears  to  increase  linearly 
with  precompression.   In  addition,  the  influence  of  precom- 
pression  upon  the  modulus  of  rupture  is  significant  with 
high  strength  mortar.   The  following  empirical  relations 
are  proposed: 


f.   =  1.250  f   +  27  (3.3) 

bo  a 


for  1:1:6  mortar,  and, 
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fbo  =  1'725  fa  +  42  (3-4) 

for  1:1:4.5  mortar. 

Here,  also,  the  specifications  are  conservative  in 
limiting  the  allowable  stress  to  23  psi. 

3.5   Combined  Compression,  Bending,  and  Shear:   Various 
test  arrangements  were  devised  to  produce  several  combina- 
tions of  bending  and  shear  under  a  constant  precompression 
of  300  psi.   These  testing  arrangements  are  shown  schemat- 
ically in  Figure  3.7,  and  the  results  are  summarized  in 
Table  3.5.   Knowing  the  mean  value  of  modulus  of  rupture 
at  a  precompression  of  300  psi  (f,  ) ,  also  the  shear  strength 
(f   )  without  flexure,  and  the  combined  bending  and  shearing 

stresses  (f,  and  f  )  ,  the  values  of  f./f,   and  f  /f    can  be 
b      s  b'  bo      s'  so 

calculated  and  plotted  to  provide  the  unit  interaction  curve 
shown  in  Figure  3.8.   It  is  anticipated  that  this  curve  is 
one  from  a  family  of  curves  that  will  be  generated  by  vary- 
ing the  precompression  in  an  expanded  program  of  testing. 
Thus,  the  flexural  and  shearing  stresses  in  a  mortar  joint 
under  combined  loading  can  be  assumed  to  follow  the  equation: 


so;     *■  boJ 


CHAPTER  4 
MODEL  TESTS  OF  CONCRETE  MASONRY  PIERS 

4-l   Selection  of  Model:   The  problem  under  investigation 
consists  of  determining  the  stress  distribution  in  a  con- 
crete masonry  shear  wall,  particularly  in  the  region  be- 
tween openings  called  piers.   Even  if  models  were  used, 
it  is  not  economical  to  fabricate  a  model  for  a  multiple 
story  wall  with  a  large  number  of  openings.   Schneider  (46) 
conducted  a  photoelastic  analysis  on  plastic  models  with 
configurations  as  shown  in  Figure  4.1.   The  loads  were 
applied  either  diagonally  or  horizontally  for  all  config- 
urations except  in  the  case  of  multiple  pier  models  where 
the  load  was  applied  in  a  manner  simulating  the  way  in 
which  a  roof  diaphragm  would  load  a  shear  wall.   A  simi- 
larity in  fringe  patterns  was  found  to  exist  between  the 
model  walls  containing  several  piers  and  a  diagonally 
loaded  single  pier.   Accordingly,  the  test  specimen  was 
chosen  as  described  in  subsection  1.3.9  of  Chapter  1. 

In  this  investigation,  a  linear  plane  stress  finite 
element  analysis  is  adopted  in  lieu  of  a  photoelastic 
analysis  to  verify  the  above  findings  and  also  to  check 
the  effect  of  pipe  columns  that  have  been  used  to  maintain 
the  geometry  of  the  openings.   The  configurations  chosen 
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Fig.  4.1.   Photoelastic  model  configurations 
chosen  by  Schneider. 
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for  analysis  and  the  finite  element  idealization  are 
shown  in  Figures  4.2  through  4.5.   Figure  4.4  shows  the 
pier  in  a  story  with  openings  and  walls  on  either  side  of 
it,  but  loaded  diagonally;  Figure  4.5  shows  the  configu- 
ration adopted  by  Schneider.   In  this  case,  the  pipe 
columns  are  idealized  as  an  equivalent  rectangular  element. 
The  equivalent  modulus  of  elasticity  of  the  element  is 
calculated  by  equating  the  axial  rigidities  of  the  actual 
pipe  column  with  that  of  the  idealized  element.   In  a 
photoelastic  analysis,  the  fringe  patterns  represent  the 
lines  joining  the  points  of  constant  maximum  shear  stress, 
also  called  isoshear  lines.   In  the  analysis,  maximum 
shear  stresses  are  computed  at  centroids  of  all  elements. 
The  contour  lines  for  these  values  represent  the  isoshear 
lines.   The  isoshear  lines  are  plotted  for  all  the  piers 
investigated  and  these  are  shown  in  Figures  4.6  through 
4.11.   It  is  evident  that  the  fringe  patterns  in  all  the 
above  cases  are  very  similar  and  this  confirms  the  method 
of  loading  adopted  in  a  single  pier  to  simulate  the  actual 
loading  conditions  a  pier  in  a  wall  is  subjected  to. 

To  check  the  effect  of  substituting  pipe  columns  for 
end  piers  (Figure  4.4),  the  modulus  of  elasticity  of  the 
pipe  column  elements  was  arbitrarily  increased  by  100%. 
However,  the  corresponding  changes  in  the  maximum  stresses 
in  the  pier  were  less  than  1%  of  the  previous  values. 
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Fig.  4.2.   Finite  element  idealization  for  a 
shear  wall  with  openings. 
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Fig.  4.3.   Finite  element  idealization  for  a 
multiple  pier  shear  wall. 
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Fig.  4.8.   Isoshear  lines  for  the  pier  shaded  as  shown  above 


80 


Fig.  4.9.   Isoshear  lines  for  the  pier  shaded  as  shown  abo\ 


"ig.  a.  10.   Tsoshear  lines  for  the  pier  shaded  as  shown  above 


82 


Fig.  4.11.   Isoshear  lines  tor  the  pier  shaded  as  shown  above. 
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On  the  basis  of  the  preceding  findings,  the  model 
configuration  shown  in  Figure  4.12  was  adopted.   The 
present  investigation  is  restricted  to  a  pier  having  a 
shear-span/depth  ratio  of  0.5. 

4.2   Selection  of  Model  Materials:   The  materials  used  in 
the  construction  of  a  model  structure  depend  upon  the 
type  of  test  to  be  made  on  the  model.   Inelastic  models 
should  be  constructed  of  materials  which  simulate  the 
physical  properties  as  closely  as  possible  in  all  respects 
(26)  (38) .   The  selection  of  a  model  material  to  simulate 
concrete  requires  the  selection  of  a  cementatious  material 
with  a  stress-strain  curve  up  to  failure  that  is  homologous 
to  that  of  concrete,  and  with  Poisson's  ratio,  ratio  of 
tensile  strength  to  compressive  strength,  and  shrinkage 
equal  to  the  counterpart  prototype  properties  (40) (9) . 
Previous  research  (40)  indicates  that  it  is  not  possible 
to  obtain  a  model  material  corresponding  to  prototype  con- 
crete by  merely  scaling  the  individual  components,  e.g., 
coarse  and  fine  aggregates,  cement,  and  the  admixtures. 
However,  as  long  as  the  overall  physical  properties  of  the 
model  material,  such  as  the  stress-strain  curve  and  the 
failure  envelope,  are  compatible  with  those  of  the  proto- 
type concrete  according  to  the  laws  of  similitude,  it  does 
not  matter  how  these  properties  are  obtained  (54) . 
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Fig.  A. 12.   Schematic  test  setup  for  finding  the 
shearing  strength  of  a  pier. 
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In  this  investigation,  model  concrete  blocks 
one-fourth  the  size  of  the  full  size  blocks  were  used. 
The  physical  and  mechanical  properties  of  the  block  have 
already  been  described  in  Chapter  2.   While  the  stress- 
strain  curves  and  strengths  of  the  model  blocks  are  com- 
parable to  the  prototype  blocks  available  in  the  market, 
their  absorption  characteristics  differ  widely  from  those 
of  the  prototype  blocks.   Because  of  their  higher  absorp- 
tion, the  water  in  the  mortar  was  absorbed  when  placed  on 
a  dry  block,  impairing  the  workability  of  the  mortar. 
Hence,  the  model  blocks  were  prewetted  before  placing 
the  mortar. 

The  types  of  mortar  used  and  their  components  are 
described  in  Section  2.6  of  Chapter  2.   To  cope  with  the 
increased  absorption  of  the  blocks  and  reduced  thickness 
of  joints  in  the  model  specimens,  an  increase  in  the 
initial  flow  of  mortar  was  needed.   This  resulted  in  re- 
duced strength  of  mortars  used  in  the  model  specimens 
compared  to  the  strengths  of  corresponding  types  of  mortars 
used  in  the  investigation  described  in  Chapter  3.   Since  it 
is  the  amount  of  aggregate  and  not  the  gradation  that  has 
the  greatest  effect  on  the  mechanical  properties  of  mortar 
(40) ,  no  attempt  was  made  to  model  the  mortar  sand  on  the 
basis  of  the  grading  curve. 

The  details  of  the  grout  mixes  used  and  their 
properties  are  given  in  Chapter  2.   The  grout  strength  and 
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its  fluidity  without  separation  of  grout  components  were 
the  factors  governing  the  choice  of  the  mixes  adopted. 

The  model  reinforcing  steel  exhibited  a  higher  yield 
strength  than  the  normal  grades  of  steel  used  in  practice. 
Since  this  investigation  is  mainly  on  nonreinforced  walls 
and  the  vertical  reinforcing  in  end  cells  do  not  contribute 
appreciably  to  the  strength  of  the  pier,  this  deviation  is 
not  considered  important.   The  depressions  in  the  bar  simu- 
late to  a  certain  extent  a  deformed  bar.   However,  no 
attempt  was  made  to  simulate  bond  characteristics.   Work- 
manship is  an  important  factor  in  masonry  which  cannot  be 
modelled. 

4.3   Fabrication  of  Model :   As  the  model  concrete  blocks 
are  one-fourth  full  size,  the  individual  pier  was  modelled 
to  represent  a  story  height  of  81  0".   Because  of  the 
limitations  of  available  space  in  the  testing  machine 
which  was  suitable  for  testing  models,  the  height  of  the 
model  was  restricted  to  2'  0"  representing  a  story  height 
of  8'  0".   Since  the  region  of  interest  is  the  pier,  it 
was  decided  to  use  cement  mortar  (1:3  by  weight)  for  the 
spandrel  portion  to  economize  the  use  of  blocks.   The 
thickness  of  mortar  joints  varied  between  0.125"  to  0.15". 
The  dimensions  of  the  model  are  shown  in  Figure  4.13. 

The  blocks  were  soaked  in  water  for  a  minute  and  then 
wiped  clean  with  paper  towels  and  allowed  to  dry  for  about 
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ten  minutes.   The  mortar  was  mixed  to  an  initial  flow  of 
about  120%.   Mortar  was  applied  to  the  web  as  well  as  the 
face  shells  of  the  block.   Also,  the  vertical  joints 
extended  to  the  full  thickness  of  the  blocks.   After  each 
course  was  laid,  the  mortar  protruding  in  the  cells  was 
removed  so  as  to  facilitate  proper  grouting  later.   After 
the  six  courses  of  blocks  were  laid,  the  pier  was  allowed 
to  cure  in  air  for  twenty- four  hours.   The  pier  was  placed 
on  a  wooden  plank  and  the  vertical  reinforcing  bars  were 
inserted  through  holes  corresponding  to  the  centers  of  the 
end  cells  in  the  pier,  and  projecting  6"  at  top  and  bottom 
of  the  pier.   The  grout  was  then  poured  in  the  cells  and 
well  tamped.   In  the  case  of  hollow  piers,  only  the  end 
cells  were  grouted.   After  being  allowed  to  cure  in  air 
for  twenty- four  hours,  this  assembly  was  placed  horizon- 
tally in  a  wooden  mold  to  form  the  2"  x  6"  x  24"  spandrels. 
The  spandrel  beam  was  reinforced  as  shown  in  Figure  4.13. 
Steel  angles  1-1/4"  x  1-1/4"  x  1/4"  two  inches  long  were 
placed  at  the  loaded  corners  to  avoid  local  crushing  of 
the  specimen.   The  spandrel  concrete  consisting  of  one 
part  of  Portland  Cement  to  three  parts  of  sand  by  weight 
with  a  water-cement  ratio  of  0.55  was  poured.   Moist  cloths 
were  kept  over  the  spandrel  surface  to  avoid  shrinkage 
cracks.   The  formwork  was  removed  after  three  days  and  the 
specimen  allowed  to  cure  in  air  until  testing. 
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Fig.  4.13.   Details  of  pier  model, 
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In  the  case  of  hollow  piers,  holes  were  left  in  the 
spandrel  beam  for  insertion  later  of  prestressing  rods 
(Figure  4.12).   This  was  accomplished  by  inserting  steel 
tubes  wound  with  greased  paper  before  pouring  the  mortar. 

4.4  Test  Setup:   The  test  setup  adopted  is  shown  in 
Figures  4.12  and  4.14.   The  pipe  columns  were  used  to 
maintain  the  geometry  of  the  openings.   They  consisted  of 
two  steel  tubes  1.0"  in  outside  diameter  and  0.05"  thick. 
Two  strain  gages  were  mounted  diametrically  opposite  to 
each  other  at  the  center  of  each  pipe  column  and  connected 
in  series  to  measure  the  load  taken  by  the  column.   The 
gages  were  previously  calibrated  for  the  strains  corre- 
sponding to  the  axial  load  on  the  column.   The  deflection 
at  the  corner  of  the  spandrel  is  measured  in  a  direction 
parallel  to  the  soffit  of  the  spandrel  beam  as  shown  in 
Figure  4.12.   This  would  correspond  to  the  horizontal  de- 
flection of  the  pier  had  the  pier  been  placed  vertically 
and  the  load  applied  diagonally. 

In  the  case  of  hollow  piers  subject  to  external 
precompression,  the  compressive  stress  was  applied  by 
tightening  two  steel  rods,  0.44  5"  in  diameter,  threaded 
at  the  ends.   Two  bearing  plates,  1/2"  x  2"  x  12",  were 
used  to  distribute  the  load  uniformly  over  the  cross- 
section  of  the  pier.   Two  strain  gages  were  mounted 
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Fig. 4. 14.    Model    test   on  pier 
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diametrically  opposite  to  each  other  on  each  rod  and 
connected  in  series  to  measure  the  axial  force  in  the  rod. 

4.5  Test  Parameters:   In  earlier  studies  by  Schneider, 
the  main  parameters  considered  were  the  shear-span  to 
depth  ratio,  the  effect  of  reinforcement  and  constant 
external  precompression,  on  different  pier  configurations. 
In  this  investigation,  only  one  pier  configuration  was 
considered.   This  study  was  designed  to: 

1)  find  the  influence  of  different  types  of  mortar  and 
grout  on  the  shear  strength  of  a  grouted  pier,  and 

2)  ascertain  the  effect  of  different  magnitudes  of 
precompression  on  the  strength  of  nongrouted  piers. 

The  mechanical  properties  of  mortar  and  grout  used 
are  described  in  Chapter  2. 

4.6  Tests  on  Grouted  Piers:   Three  different  types  of 
grout  and  two  types  of  mortars  were  used  in  this  test 
series.   In  all  the  cases,  the  failure  occurred  by  diagonal 
tensile  cracking  along  the  vertical  diagonal,  the  first 
crack  occurring  near  the  center  of  the  pier  at  about  85% 

of  the  ultimate  load.   The  diagonal  crack  then  progressed 
quickly  towards  the  corners.   In  no  case  did  complete 
collapse  occur  even  after  large  cracks  had  opened  up  com- 
pletely through  the  pier.   The  failure  patterns  of  piers 
with  strong  mortar  joints  and  different  types  of  grout 
are  shown  in  Figures  4.15  through  4.17.   The  failure 
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Fig.  4.15.   Failure  pattern  of  PR  4 

strong  mortar,  weak  grout. 


Fig.  4.16.   Failure  pattern  of  PR  3 

strong  mortar,  medium  grout, 
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Fig.  4.17.   Failure  pattern  of  PR  5 

strong  mortar,  strong  grout. 


Fig.  4.18.  Failure  pattern  of  PR  9 
weak  mortar,  weak  grout, 
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Fig.  4.19.  Failure  pattern  of  PR  10 

weak  mortar,  medium  grout. 


Fig.  4.20.  Failure  pattern  of  PR  7 

weak  mortar,  strong  grout. 
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Fig.  4.21.   Typical  load-deflection  curves  for  grouted  piers, 
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1 

Fig.  4.29.   Failure  pattern  of  hollow  pier  PR  6 
(weak  mortar) . 
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patterns  of  piers  with  weak  mortar  joints  are  shown  in 
Figures  4.18  through  4.20.   Typical  load  deflection  curves 
are  shown  in  Figure  4.21.   For  each  test,  two  corresponding 
three-block  grouted  prisms,  three  mortar  cubes,  three 
grout  cubes,  and  a  model  concrete  block  were  tested  in 
axial  compression. 

4.7   Tests  on  Nongrouted  Piers:   In  order  to  study  the 
effect  of  external  precompression  on  nongrouted  piers  with 
weak  and  comparatively  strong  mortar  joints,  piers  were 
tested  at  precompressions  of  0,  50,  125,  and  200  psi.   It 
was  found  that  in  the  case  of  piers  with  weak  mortar  joints, 
the  failure  was  always  initiated  by  debonding  of  mortar 
joints  irrespective  of  the  value  of  precompression  imposed. 
However,  the  piers  with  stronger  mortar  joints  failed  by 
diagonal  tensile  cracking  along  the  vertical  diagonal  at 
precompressions  above  50  psi.   In  all  the  cases  where 
external  precompression  was  applied,  the  joints  near  the 
center  of  the  pier  failed  first  and  the  crack  extended 
towards  the  corners.   The  failure  patterns  are  shown  in 
Figures  4.22  through  4.27.   Figures  4.28  and  4.29  show 
the  mode  of  failure  of  a  pier  with  strong  and  weak  mortar 
joints,  respectively,  when  no  external  precompression  was 
applied.   In  no  case  was  a  complete  collapse  of  the  whole 
specimen  observed  after  attainment  of  the  maximum  load. 

The  tests  results  are  presented  in  the  next  chapter. 


CHAPTER  5 
ANALYSIS  OF  TEST  RESULTS 

5.1   Grouted  Piers:   Table  5.1  presents  the  results  of  the 
tests  on  grouted  piers.   The  ultimate  shear  strength  of  the 
pier  was  calculated  as  the  average  shear  stress  over  the 
gross  area  of  the  pier  as: 


fs  =  H/tD  (5.1) 


where  H  equals  the  horizontal  component  of  the  ultimate 
load  in  lbs  (P  //?) ,  t  equals  the  thickness  of  wall  in 
inches  (1.9"),  D  equals  the  overall  depth  of  the  pier  in 
inches  (12"),  and  f   equals  the  average  shear  stress  over 
the  cross-section  of  the  pier  (psi). 

Figure  5.1  shows  the  plot  between  the  square  root  of 
the  prism  strength  (f  ' )  and  the  average  shear  stress  on 
the  pier.   Using  the  method  of  least  squares,  the  best-fit 
curve  for  the  observed  data  was  found  to  plot  as  a  straight 
line  yielding  the  following  relation: 


f   =  3.8  /f~7  (5.2) 

s         m 
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WEAK    MORTAR 
STRONG    MORTAR 
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Fig.  5.1 


Relation  between  average  shear  stress  and  square  root  of 
prism  strength  for  grouted  piers. 


105 


Figure  5.2  shows  the  plot  between  f  '  and  f  which  yields 

m       s       J 

the  relation: 


fs  =  20  +  0.0829  fm'  (5.3) 


It  is  not  possible  to  attach  any  physical  significance 
between  the  failure  of  a  grouted  prism  in  axial  compression 
and  a  grouted  pier  in  diagonal  tension.   However,  so  long  as 
the  mortar  and  grout  are  of  strengths  used  in  this  investi- 
gation, Equation  (5.2)  could  be  used  as  an  empirical  equa- 
tion to  assess  the  diagonal  tensile  strength  of  grouted 
masonry  piers. 

It  can  also  be  noted  that  the  strengths  of  the  piers 
with  weak  and  strong  mortar  joints,  but  with  the  same 
strength  of  grout,  are  not  appreciably  different. 

5.2  Nongrouted  Piers:   The  ultimate  shear  stress  in  the 
case  of  a  nongrouted  pier  is  calculated  on  the  basis  of  the 
net  area  of  the  cross-section  shown  shaded  in  Figure  5.3. 
The  applied  precompression  was  also  based  on  this  area  only. 
Since  the  load  is  applied  diagonally,  the  vertical  component 
of  the  load  would  contribute  an  internal  compressive  stress 
in  addition  to  the  external  precompression  applied.   The 
internal  precompressive  force  is  the  vertical  component  of 
the  load  less  the  load  taken  by  the  pipe  columns.   It  is 
assumed  that  both  the  compressive  stresses  are  uniformly 
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distributed  over  the  cross-section  of  the  pier.   The 
results  of  the  tests  on  nongrouted  piers  are  shown  in 
Table  5.2.   Figure  5.4  shows  the  plot  between  the  total 
precompression  and  the  average  shear  stress  in  the  pier. 
It  is  evident  that  precompression  definitely  aids  in 
increasing  the  shear  capacity  of  a  nongrouted  pier.   Also, 
there  was  not  a  marked  increase  in  the  shear  capacity  of 
the  piers  when  stronger  mortar  joints  were  used.   However, 
this  observation  is  limited  to  the  strengths  of  mortars 
used  in  this  investigation. 

5. 3   Equations  for  Predicting  the  Diagonal  Tensile  Strength 
of  Masonry: 

5.3.1   Grouted  Piers:   on  the  basis  of  this  investigation, 
it  is  concluded  that  Equation  (5.2)  could  be  used  for  pre- 
dicting the  diagonal  tensile  strength  of  grouted  masonry 
for  the  ranges  of  strengths  of  mortar  and  grout  and  for  the 
pier  configuration  considered  in  this  investigation.   Blume 
(6)  has  suggested  that  diagonal  tensile  strengths  should 
always  be  ascertained  by  testing  a  corresponding  square 
panel  under  diagonal  compression.   In  a  model  investigation 
such  as  this  one,  it  is  difficult  to  define  a  suitable  size 
for  the  companion  panel.   Also,  results  on  4-ft-square 
panels  tested  by  Schneider  indicated  that  the  ultimate 
shear  strength  of  the  masonry  was  nearly  the  same  irrespec- 
tive of  the  reinforcing  present  in  these  panels  and  it  was 
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not  the  same  as  the  diagonal  tensile  stress  produced  by 
combined  flexural  tension  and  shear  stresses  in  a  pier 
contained  within  a  wall  (46) .   In  the  light  of  the  above 
remarks,  it  is  suggested  that  Equation  (5.2)  is  a  suitable 
means  for  ascertaining  the  strength  of  grouted  piers  with 
shear  span  to  depth  ratio  (a/D)  around  0.5.   Since  the 
shear  strength  varies  with  a/D  ratios,  more  tests  are 
necessary  to  arrive  at  a  suitable  value  corresponding  to 
other  pier  configurations. 

5.3.2   Nongrouted  Piers:   a  pier  subjected  to  diagonal 
compression  is  idealized  as  shown  in  Figure  5.5.   Applying 
the  simple  beam  theory  and  considering  YY  as  the  neutral 
plane,  the  bending  stress  at  any  section  XX  is: 


fb  =  (Hy  -  XjVJ/S 


y.-A 


where  S    is  the  elastic  section  modulus  corresponding  to 
the  extreme  fiber  computed  on  the  basis  of  the  cross-section 
shown  in  Figure  5.3,  and  H  and  V  are  the  horizontal  and 
vertical  components  of  the  diagonal  load  'P'.   The  total 

v-v 

precompression  on  the  cross-section  is  taken  as  ( — ^ — 
+  fpr)  where  A  is  the  area  of  cross-section  of  the  pier  as 
shown  in  Figure  5.3,  f    is  the  extent  of  external  precom- 
pression applied,  and  VpT  is  the  load  taken  by  the  pipe 
columns.   Using  Equation  (3.3) ,  the  maximum  allowable 
flexural  tension  at  this  level  of  precompression  can  be 
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Fig.  5.6.   Test  results  of  nongrouted  pier  speciemens . 
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found  (fbQ) •   The  average  shear  stress  on  the  cross-section 
is  calculated  as: 


fs  =  H/A 


and  the  shear  strength  of  the  cross-section  at  this  level 

of  precompression  is  calculated  from  Equation  (3.1)  as  f 

so 

It  was  verified  whether  the  combined  shearing  and  bending 
stresses  on  the  mortar  joint  satisfy  the  interaction  equa- 
tion proposed  in  Chapter  3,  viz., 

v  soJ     v  bo-' 

The  above  results  are  presented  in  Table  5.3.   It  can 

be  seen  that  the  scatter  is  excessive  in  the  case  of  piers 

with  weak  mortar  joints.   However,  other  pier  specimens 

seem  to  follow  the  interaction  equation  quite  closely.   The 

values  of  f  /f    and  f,  /f,  „  are  plotted  on  the  interaction 
s   so      b  bo     c 

diagram  to  show  the  deviations  graphically,  as  shown  in 
Figure  5.6.   The  observed  mode  of  failure  verifies  the  above 
equation  only  for  the  piers  without  precompression  where  the 
debonding  occurred  first  near  point  '3'  in  Figure  5.5.   In 
all  other  cases,  the  crack  originated  at  or  near  the  center 
of  the  panel.   This  suggests  that  higher  precompression 
resulted  in  correspondingly  higher  flexural  tensile  strength 
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at  point  '3',  thus  suppressing  debonding  of  the  mortal"  joint 
at  that  point.   Hence,  it  is  logical  to  assume  that  the 
critical  point  under  higher  precompression  for  a  pier  config- 
uration considered  in  this  investigation  would  be  the  center 
of  the  pier  wherein  the  bending  stresses  are  absent.   There- 
fore, it  can  be  concluded  that  failure  by  debonding  would 
occur  in  those  cases  wherein  the  shear  stress  exceeds  the 
shear  strength  of  the  joint  given  by  Equations  (3.1)  and 
(3.2),  depending  on  the  type  of  mortar  used.   This  would 

suggest  a  value  of  f  /f    =  1  in  column  (9)  of  Table  5.3. 
^  s   so 

As  can  be  seen,  the  values  of  f  /f  Q  range  from  0.7  to  0.916, 
suggesting  the  variability  involved.   It  can  also  be  noted 
that  Equations  (3.1)  and  (3.2)  consistently  overestimate  the 
capacity  of  the  pier.   Hence,  it  is  suggested  that  the 
shear  strength  obtained  using  the  above  equations  be  suitably 
factored.   A  multiplier  of  0.7  is  recommended.   The  adjusted 
values  using  a  factor  of  0.7  are  shown  in  column  (15)  of 
Table  5.3. 

The  above  procedure  is  recommended  for  piers  with  an 
a/D  ratio  less  than  or  equal  to  0.5.   For  higher  a/D  ratios, 
the  bending  stresses  are  more  predominant  and  the  interaction 
Equation  (3.5)  seems  to  be  more  applicable. 

Since  the  debonding  failure  mode  could  be  suppressed 
completely  at  higher  levels  of  precompression,  it  can  be 
conceived  that  failure  might  be  initiated  by  exceeding  the 
tensile  strength  of  the  mortar  or  block.   The  maximum 
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principal  stress  occurs  at  the  center  of  the  pier.   The 
state  of  stress  at  this  point  is  assumed  to  be  as  shown  in 
Figure  5.8.   The  maximum  principal  stress  is  given  by: 


f    =  -  I(V  +  f   )  +  /u725(f  +  f   )  2  +  (|)  2 
max      2VA     pr'         VA     pr'     VA' 


Tensile  cracking  is  assumed  to  occur  if  f    exceeds  the 
3  max 

tensile  strength  of  the  block  or  mortar,  whichever  is  less: 


f     >  f_  (5.4) 

max  —   t 


The  values  of  f    /f.  are  given  in  column  (13)  in  Table  5.3. 
max   t 

The  results  indicate  that  in  no  case  did  failure  occur  due 
to  tensile  cracking. 

Another  failure  criterion  that  could  be  conceived  is 
to  equate  the  maximum  principal  stress  to  the  diagonal  ten- 
sile strength  of  a  square  panel.   A  square  panel  8"  x  8" 
was  built  for  this  purpose  with  the  same  mortar  and  blocks 
used  for  building  the  corresponding  pier.   This  was  tested 
under  diagonal  compression  as  shown  schematically  in 
Figure  5.7.   For  both  mortar  types,  the  failure  was  initi- 
a1  3     debonding  of  the  mortar  joints  as  shown  in  Figure 

panel  witl  v  -aker  mortar  joints  failed  at  3  50  lbs, 
while  the  one  with  stronger  joints  failed  at  600  lbs. 
Based  on  the  net  area  of  cross-section  of  the  panel,  the 
average  shear  stress  as  a  measure  of  the  diagonal  tension 
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NET   AREA   =  8.6  SQ0 IN. 


Fig.  5.7.   Failure  pattern  of  hollow  square  panel  subject  to 
diagonal  compression. 
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V/A  +  L 


Fig.  5.8.   Assumed  state  of  stress  at  the  center  of  pier, 
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for  both  cases  is  ,  respectively,  29  psi  and  49  psi.   The 

values  of  f   /f ,.  are  given  in  column  (14)  in  Table  5.3, 
max  at 

where  f ,.  is  the  value  of  the  shear  stress  obtained  in 
at 

the  diagonal  tension  test.  As  the  results  indicate,  this 
criterion  seems  to  be  valid  only  for  the  case  without  any 
external  precompression. 

The  equations  proposed  by  Sinha  and  Hendry  (50)  for 
predicting  the  load  capacity  of  brickwork  have  already  been 
presented  in  Chapter  1.   They  proposed  that  within  certain 
values  of  precompression  the  failure  was  initiated  by  the 
maximum  stress  exceeding  the  tensile  strength  of  the 
brickwork.   In  the  case  of  brickwork,  the  tensile  strength 
was  determined  by  a  diametrical  compression  test  on  a  circu- 
lar disc  of  masonry.   The  equations  proposed  by  Sinha  and 
Hendry  (50)  were  applied  for  the  nongrouted  piers  assuming 
the  diagonal  tensile  strength  obtained  from  the  square 
panel  for  f   in  Equation  (1.10)  and  Equation  (1.11). 
Imaginary  values  were  obtained  for  the  lower  limit  of  pre- 
compression.  The  upper  limits  of  precompression  were 
found  to  be  85  psi  and  158  psi  for  piers  with  weak  and 
strong  mortar,  respectively.   This  implies  that  in  the  case 
of  piers  with  weak  mortar  joints,  the  failure  will  result 
in  compression  typically  by  vertical  splitting  for  precom- 
pressive  stresses  above  8!i  psi.   In  this  study,  the 
failure  patterns  suggest  a  failure  initiated  by  debonding 
of  mortar  joints.   The  difference  may  be  attributed  to  the 
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size  of  the  square  panel  considered  for  determining  the 
diagonal  tensile  strength  and  in  differences  between  brick 
and  concrete  block  masonry.   However,  more  tests  are  nec- 
essary to  arrive  at  a  definite  conclusion  in  this  regard. 

In  the  light  of  the  preceding  remarks,  the  following 
equations  are  proposed  for  ascertaining  the  shear  strength 
of  a  nongrouted  masonry  pier: 

Case  (i) :   The  following  equations  are  considered  to 
be  applicable  for  pier  configurations  for  low  levels  of 
precompression;  but  for  piers  with  shear  span  to  depth  ratio 
>  0.5,  these  equations  may  be  considered  to  be  applicable 
for  all  levels  of  precompression: 

*  so-*     *■  bo-' 

f    =  f  i. 

max    t 

The  lesser  of  the  loads  predicted  by  the  above  equations 
is  to  be  taken  as  the  ultimate  load  on  the  pier. 

Case  (ii) :   For  piers  with  a/D  ratio  less  than  0.5, 
the  governing  equations  are: 

f   =  0.7(f   ) 
s        so 

f    =  f . 

max     t 

The  minimum  load  obtained  from  the  above  equations  is  to 
be  taken  as  the  ultimate  load  on  the  pier. 


CHAPTER  6 
NONLINEAR  FINITE  ELEMENT  ANALYSIS 

6.1   General  Remarks :   The  finite  element  method  is  now 
widely  accepted  as  a  method  of  stress  analysis.   "Progress 
in  the  method  has  been  made  on  three  fronts,  all  of  which 
contribute  to  the  strength  and  flexibility  of  the  method. 
First  of  all,  the  relation  of  the  finite  element  method  to 
previous  well-established  methods  in  continuum  mechanics 
has  given  it  a  firm  foundation.   Secondly,  the  search  for 
and  development  of  the  many  consistent  elements  has  given 
it  a  wide  area  of  application.   Finally,  extension  of  the 
method  to  the  study  of  nonlinear  behavior  in  both  material 
and  geometric  nonlinearities  has  resulted  in  more  realistic 
models  and  design  methods"  (37,  p.  35). 

The  finite  element  method  has  been  described  extensively 
in  the  literature.   A  comprehensive  discussion  of  the  theory 
and  application  of  the  method  to  two-dimensional  and  three  - 
dimensional  systems  is  given  in  the  book  by  Zienkiewicz  (55) . 
An  introductory  treatment  of  the  method  is  given  in  the  book 
by  Desai  and  Abel  (15) . 
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6.2  Finite  Element  Linear  Analysis;   The  following  steps 
outline  the  basic  procedure  of  a  Finite  Element  Linear 
Analysis : 

1.  Discretization  of  the  continuum:   The  basis  of  the 
finite  element  method  is  the  representation  of  a  body  or  a 
structure  by  an  assemblage  of  subdivisions  called  finite 
elements.   These  elements  are  interconnected  at  joints,  which 
are  called  nodes,  and  also  along  their  interfaces. 

2.  Selection  of  displacement  model:   Simple  functions 
called  interpolation  functions  are  chosen  to  approximate  the 
distribution  or  variation  of  the  actual  displacements  (or 
other  field  variables)  over  each  finite  element.   Thus,  the 
element  displacements  'u(x,y,z)'  are  expressed  in  terms  of 

a  finite  number  of  nodal  displacements  (<5)  through  assumed 
interpolation  functions  'N': 

(u(x,y,z)}  =  {N(x,y,z)}  {6}  (6.1) 

It  is  generally  not  possible  to  select  a  displacement  func- 
tion that  can  represent  exactly  the  actual  variation  of  dis- 
placement in  the  element.   Hence,  the  basic  approximation 
of  the  finite  element  method  is  introduced  at  this  stage. 

3.  Strains:  By  differentiating  equation  (6.1),  ele- 
ment strains  {e}  may  be  stated  in  terms  of  the  nodal  point 
displacements  {&}: 

is]    =    {B}  {&}  (6.2) 
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4.   Stresses:   Stresses  are  related  to  strains  through 
a  constitutive  law.   Usually,  a  linear  relation  is  assumed. 


{a}   =  [D]  ({el  -  (ec})  +  {%} 


(6.3) 


where  {e  }  are  the  initial  strains  and  {oQ}    the  residual 
stresses. 

5.  Using  the  principle  of  virtual  displacements,  the 
element  stiffness  matrix  [k]  relating  element  nodal  forces 
{F}  to  nodal  displacements  {6}  is  formed: 

{F}  =  [k]  {6}  (6-4) 

where, 

[k]  =  f  (B}T  [D]  {B}  d(vol)  (6.5) 

6.  The  structure  stiffness  matrix  or  global  stiffness 
matrix  [K]  relating  external  nodal  point  forces  {R}  to  the 
nodal  displacements  {6}  is  assembled  by  using  the  Direct 
Stiffness  Method  where  it  is  possible  to  add  directly  the 
individual  stiffnesses  [k]  and  loads  to  locations  in  the 
overall  matrices  [K]  and  {R}  in  conformity  with  the  require- 
ments of  one-to-one  correspondence  between  the  nodes  of  the 
elements  and  those  of  the  assemblage. 
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{R}  =  [K]  {6} 


(6.6) 


7.   Equation  (6.6)  is  solved  to  determine  the  nodal 
point  displacements : 


{6}  =  [K]"1  {R} 


(6.7) 


8.   The  element  strains  and  the  element  stresses  are 
found  using  equations  (6.2)  and  (6.3),  respectively. 

6 . 3   Finite  Element  Nonlinear  Analysis: 

6.3.1   General  Remarks:   in  the  displacement  method  of 
finite  element  analysis,  nonlinearities  occur  in  two  dif- 
ferent forms.   The  first  is  material  or  physical  nonlinear- 
ity  which  results  from  nonlinear  constitutive  laws.   The 
second  is  geometric  nonlinearity  which  derives  from 
finite  changes  in  the  geometry  of  the  deforming  body. 

The  first  category  encompasses  problems  in  which  the 
stresses  are  not  linearly  proportional  to  the  strains,  but 
in  which  only  small  displacements  and  small  strains  are 
considered  and  linear  strain  displacement  relations  are 
used.   Since  displacements  in  reinforced  concrete  structures 
are  inherently  small,  geometric  nonlinearities  generally 
do  not  occur  and  attention  can  be  restricted  to  the  problems 
of  nonlinear  material  properties  (47).  These  include: 
1)  cracking  of  the  concrete;  2)  nonlinear  stress-strain 


126 


relations  for  concrete,  steel,  bond  and  aggregate 
interlock;  and  3)  time  dependent  effects  such  as  creep, 
shrinkage,  temperature  and  load  history. 

6.3.2   General  Physical  Approach:   for  arriving  at  the 
solution  of  a  linear  elastic  problem,  the  set  of 
equations, 

[K]  {6}  -  {R}  -  0  (6.8) 

were  solved  in  which  the  vector  {R}  lists  all  the  forces 
due  to  external  loads,  initial  stresses  and  strains,  etc. 

In  the  above,  a  linear  constitutive  law  (Eq.  6.3)  and 
a  linear  strain  displacement  relation  (Eq.  6.2)  were 
assumed.   Continuity  of  displacements  was  required  and  the 
equilibrium  conditions  were  satisfied  approximately. 

In  a  nonlinear  constitutive  problem,  displacement 
continuity  and  equilibrium  still  have  to  be  satisfied.   Thus, 
the  only  relationship  which  has  to  be  replaced  is  Equation 
(6.3).   This  new  relation,  in  general,  will  be  given  by  some 
relation, 

f  {{a},    {e})  =  0  (6.9) 

If  a  solution  of  Equation  (6.8)  can  be  achieved  in  which, 
by  adjustment  of  one  or  more  of  the  parameters  [D] ,  {tQ}>    or 
{<:•    )  in  Equation  (6.3),  this  and  Equation  (6.9)  are  made  to 
yield  the  same  stress  values,  then  a  solution  is  found.   An 
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iterative  procedure  becomes  necessary.   If  the  iteration  is 
conducted  by  adjustment  of  the  [D]  matrix,  the  process  is 
known  as  one  of  variable  stiffness.   If  ieQ}   or  {oq}  is 
adjusted,  so-called  initial  strain  or  initial  stress 
approaches  are  obtained. 

In  many  real  situations,  relations  of  the  type  of 
Equation  (6.9)  cannot  be  written  in  terms  of  total  strains 
and  stresses;  however,  relations  could  be  established  be- 
tween the  increments  {Aa}  and  Ue}.   In  such  cases,  the 
same  process  can  be  applied  for  an  increment  of  load. 
Such  incremental  processes  can  be  combined  with  any  of  the 
previous  methods. 

6-4   Variable  Stiffness  Methods:   The  procedure  for  linear 
analysis  described  in  Section  6.2  can  be  utilized  in  con- 
junction with  one  of  several  iterative  schemes  to  account 
for  nonlinear  behavior.   The  simple  truss  structure  of 
Figure  6.1  can  be  utilized  to  explain  the  three  common 
iterative  solution  techniques  illustrated  in  Figure  6.2  (47) 
These  are  designated  as:   a)  initial  stiffness  method, 
hs      ecant  stiffness  method,  and  c)  tangent  stiffness  method. 

<-roSs-strain  relation  for  any  of  the  elements  is 

r.o  be  nonl:ne'-  in  compression  and  to  be  linear  in 
tension.  The  nonlinear  analysis  is  to  be  made  for  a  given 
load  R.   The  following  steps  are  performed: 

1.   For  the  segment  0-1  in  Figure  6.2,  the  initial 
stiffness  0-1  on  the  o-e  curve  is  used  in  all  the  three 
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Fig.    6.1.      Effect   of  nonlinear  material  properties. 
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Fig.    6.2.      Iterative  solution  techniques   for   nonlinear 
analysis. 
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methods  to  find  the  eleme>nt  stiffness  [k]  and  the  structure 
stiffness  [K] .   Equations  (6.1)  through  (6.7)  are  used  to 
find  the  displacements  and  strains. 

2.  For  the  segment  1-2  in  all  three  methods,  the 
actual  stress  a  corresponding  to  e  from  step  1  is  found  for 
each  element.   The  unbalanced  element  stresses  Aa, „  are  com- 
puted, from  which  the  unbalanced  nodal  forces  AR.„  can  be 
obtained  by  proper  summation. 

3.  For  the  segments  2-3  and  3-4,  the  procedure  differs 
for  the  three  methods.  In  each  case,  the  structure  is  to  be 
analyzed  for  the  unbalanced  load  AR,  2« 

a)  the  initial  stiffness  from  the  a-e  curve  is 
used  again  and  steps  1  and  2  above  are  repeated. 

b)  the  secant  stiffness,  0-2-3,  from  the  a-e 
curve  is  used  and  steps  1  and  2  above  are  repeated. 

c)  the  tangent  stiffness  at  point  2  on  the  a-e 
curve  is  used  and  steps  1  and  2  above  are  repeated. 

In  all  the  three  methods,  if  a  tensile  failure  occurs 
in  an  element,  its  stiffness  [k]  is  disregarded. 

4.  Steps  1,  2,  and  3  above  are  repeated  until  the  un- 
balanced nodal  point  forces  AR  reach  a  predefined  acceptable 
value,  at  which  time,  final  values  of  stresses,  strains  and 
displacements  are  output. 
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It  is  evident  that  the  convergence  is  fastest  by  the 
tangent  stiffness  method  and  slowest  by  the  initial  stiff- 
ness method,   On  the  other  hand,  only  the  initial  stiffness 
method  possesses  the  distinct  advantage  that  the  structure 
stiffness  [K]  need  not  be  reformulated  for  each  iteration. 
If  direct  solution  methods  are  used  in  the  elastic  program, 
the  variable  stiffness  method  is  extremely  uneconomical. 


6.5   Initial  Stress  Method:   The  "initial  stress"  process 
as  proposed  by  Zienkiewicz  et  al.  (57)  approaches  the  solu- 
tion of  a  nonlinear  problem  as  a  series  of  approximations. 
In  the  first  place,  during  a  load  increment,  a  purely  elas- 
tic problem  is  solved  and  incremental  strains  {Ae1}  and 
incremental  stresses  {Ao  }  are  obtained.   The  nonlinearity 
implies,  however,  that  for  the  increment  of  strain  found, 
the  stress  increment  will.,  in  general,  not  be  correct.   If 
the  true  increment  of  stress  possible  for  the  given  strain 
is  fAa},  then  the  situation  can  only  be  maintained  by  a  set 
of  body  forces  equilibrating  the  "initial  stress"  system 
({AoM  -  {Ao}). 

At  the  second  stage  of  the  computation,  this  body 
force  system  can  be  removed  by  allowing  the  structure  (with 
unchanged  elastic  properties )  to  deform  further.   An  addi- 
tional  set  of  strain  and  corresponding  stress  increments  is 
caused.   Once  again,  these  are  likely  to  exceed  those  per- 
missible by  the  nonlinear  relationship,  and  the  redistribu- 
fcicr.  of  equilibrating  bo  iy  forces  has  to  be  repeated. 
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If  the  process  converger-'   then  finally  within  an  increment 
the  full  nonlinear  comp,  I  ibility  and  equilibrium  conditions 
will  be  satisfied,  just:  as  they  are  in  an  "incremental 
elasticity"  solution. 

It  can  be  noted  that  at  every  stage  of  the  iteration 
the  same  stiffness  rnatri::  is  used  and,  if  this  is  once 
inverted,  each  step  can  be  accomplished  in  a  small  fraction 
of  the  time  needed  foi  the  first  solution.   The  body  forces 
calculated  at  the  last  iteration  can  be  considered  as  an 
unbalanced  residual  force  left  on  the  structure  and  this 
serves  as  a  means  to  assess  the  error  involved. 

6'6   Initial  Strain  Method:   This  method  is  similar  to  the 
init  ial  stress  approach  and  is  useful  when  the  stress  level 
cannot  be  explicitly  determined  in  terms  of  strain,  but 
strains  in  terms  of  stress.   Thus  symbolically, 


(el  =  f  ({a}) 


(6.10) 


Now  the  equality  between  the  above  equation  and  the  elastic 
relation  (Eq.  6.3)  can  be  obtained  by  adjusting  {eq}.   Here, 
also,  Equation  (6.8)  is  solved  iteratively,  but  the  elastic 
strains  obtained  at  every  stage  are  compared  with  those 
cor respondi nq  to  the  constitutive  relation  (Eq.  6.10)  and 
the  difference  used  in  evaluating  the  residual  forces. 
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6.7   Previous  Investigations:   The  various  aspects  of 
application  of  the  finite  element  method  to  reinforced 
concrete  structures  were  reviewed  recently  by  Scordelis 
(47) .   The  response  of  a  reinforced  concrete  wall  panel 
specimen  was  studied  experimentally  and  analytically  by 
Cervenka  and  Gerstle  (10) .   Constant  strain  triangular 
finite  elements  were  used  to  model  the  concrete  and  steel 
reinforcement.   The  steel  was  assumed  to  be  uniformly 
distributed  over  the  element  in  two  orthogonal  directions, 
and  a  composite  material  constitutive  relation  was  devel- 
oped in  order  to  determine  the  element  stiffness  matrix. 
The  constitutive  relations  are  shown  in  Figures  6.3(a) 
and  6.3(b).   The  assumed  biaxial  stress  criterion  is  shown 
in  Figure  6.3(c).   The  concrete  plasticity  was  limited 
by  maximum  concrete  strain  e  .   Plasticity  of  the  uncracked 
concrete,  which  is  in  a  biaxial  state  of  stress,  was 
assumed  to  obey  the  von  Mises  yield  criterion  and  associ- 
ated flow  rule.   A  method  similar  to  that  presented  by 
Zienkiewicz  et  al.  (57)  was  used  in  the  finite  element 
analysis  to  account  for  plastic  deformations  under  biaxial 
compressive  yielding.   An  incremental  nonlinear  analysis 
was  used  which  accounted  for  tensile  cracking  and  subse- 
quent stress  redistribution.   No  account  was  taken  of  bond, 
the  effects  of  dowel  action  and  aggregate  interlocking  of 
crack  faces.   The  experimental  and  analytical  crack  patterns 
obtained  are  shown  in  Fi  jures  6.3(d)  and  (e)  . 
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(a)      Uniaxial   Streaa-Stralr 
Curve    for  Concrete 


-  ■ -«5, 

=LT 

^ASSUMED 


<b)      Uniaxial  Stress-Strain 
Curve    for  Steel 


(c)      Concrete    Failure  Criteria 


(d)    Experimental    Crack    Pattern 


(e)      Analytical    Crack    Patte 


Fig.  6.3.   Material  idealization  and  selected  results  of 
analysis  of  RC  panels  by  Cervenka  and  Gerstle. 
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Franklin  (19)  used  quadrilateral  discrete  elements 

for  concrete  and  one  dimensional  bar  elements  for  steel 

reinforcement  for  analyzing  shear  panels.   A  multilinear 

stress-strain  curve  was  used  for  concrete  in  compression 

and  also  for  reinforcement  to  approximate  the  actual 

curves  obtained  in  uniaxial  tests,  while  a  linear  curve 

was  assumed  up  to  failure  in  tension  for  concrete. 

Figure  6.4  shows  the  biaxial  failure  criterion  assumed 

for  concrete.   Tensile  cracking  was  accounted  for  and  the 

concrete  was  assumed  to  have  failed  when  it  reached  f  ' 

c 

in  compression  with  no  subsequent  plasticity.   An  incre- 
mental iterative  secant  stiffness  approach  was  used  for 
the  analysis. 

Yuzugullu  and  Schnobrich  (53)  used  quadrilateral 
elements  composed  of  four  constant  strain  triangles  to 
model  a  shear-wall  frame  system.   Special  tie  link  elements 
were  used  to  connect  the  wall  to  the  frame.   The  steel  bar 
reinforcement  mesh  used  in  the  wall  was  assumed  to  be 
smeared  out  and  uniformly  distributed  over  each  quadri- 
lateral element.   The  von  Mises  yield  criterion  was  assumed 
for  the  biaxial  compression  of  concrete.   An  elastic  per- 
fectly plastic  stress-strain  relationship  was  assumed  for 
concrete  and  steel,  with  a  limiting  compressive  strain  in 
the  case  of  concrete.   A  composite  material  property  matrix 
for  the  reinforced  concrete  model  was  obtained  as  the  sum 
of  the  material  property  matrices  of  the  concrete  and  of 
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Late  rol 
<i  t  r  e  n  s      t. 


Fig.    6.4.      Failure  envelope  adopted  by  Franklin. 
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reinforcing.   Full  bond  was  assumed  to  remain  even 
throughout  the  crack  region.   Concrete  was  assumed  to 
crack  when  principal  tensile  stress  reached  the  tensile 
strength  of  the  concrete.   It  was  assumed  that  the  cracked 
concrete  could  no  longer  carry  any  tensile  forces  perpen- 
dicular to  the  crack,  but  could  maintain  some  amount  of 
shear  stiffness  because  of  the  irregular  surface  of  the 
crack.   A  shear  retention  factor  3  was  introduced  so  that 
a  shear  force  could  be  transferred  across  the  cracked 
planes.   It  was  realized  that  the  shear  strength  along  the 
crack  should  be  a  function  of  the  crack  width  and  has 
upper  and  lower  bounds  of  one  and  zero  relative  to  the  un- 
cracked  shear  strength  capacity  (23).   However,  a  constant 
value  of  0.4  was  assumed  in  the  analysis.   Elastic  behavior 
of  a  cracked  element  was  terminated  if  compressive  yielding 
occurred  in  the  cracked  concrete  and/or  yielding  occurred 
in  the  steel  reinforcement.   Once  concrete  reached  its 
ultimate  strain,  it  was  assumed  to  crack  and  its  stiffness 
was  disregarded. 

Figures  6.5  and  6.6  illustrate  schematically  the 
method  of  solution  adopted.   Method  1  was  used  to  account 
for  plasticity  of  the  reinforcement  and  cracking  of  concrete, 
while  Method  2  was  used  for  biaxial  plasticity  of  the 
concrete. 

Suidan  and  Schnobrich  (52)  used  three-dimensional 
20-node  isoparametric  elements  to  investigate  the  cracking, 
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Fig.  6.5.   Schematic  diagram  to  illustrate  crack  propagation 
(initial  stress  method  using  variable  stiffness 
within  an  increment) . 


Fig.  6.6.   Schematic  diagram  to  illustrate  crack  propagation 
(initial  stress  method  using  constant  stiffness 
within  an  increment) . 


138 


crushing,  and  yielding  of  reinforced  concrete  beams.   The 
constitutive  relations  and  the  method  of  analysis  were  the 
same  as  in  the  previous  case.   However,  stresses  and  strains 
were  computed  at  integration  points  so  that  the  failure 
patterns  could  be  investigated  at  these  points.   In  both 
of  the  previous  two  investigations,  no  biaxial  stress 
failure  criteria  were  assumed. 

6.8   Objective  and  Scope  of  Present  Investigation:   In  this 
investigation,  an  attempt  is  being  made  to  use  the  Finite 
Element  Method  to  estimate  the  load  capacity  of  concrete 
masonry  elements.   The  analytical  model  and  the  method  of 
solution  adopted  are  described  in  Chapters  7  and  8.   The 
analysis  takes  into  account  cracking,  crushing,  and  com- 
pressive yielding  of  the  concrete  blocks,  mortar,  and 
grout,  as  well  as   the  debonding  of  mortar  joints.   The 
theoretical  results  are  compared  with  experimental  data. 
The  analysis  considers  only  short  time  loading. 


CHAPTER  7 
DESCRIPTION  OF  ANALYTICAL  MODEL 

7.1   Choice  of  Finite  Elements:   The  problem  under 
investigation  is  assumed  to  be  approximated  by  plane  stress 
conditions.   A  number  of  finite  elements  have  behavior 
characteristics  that  are  adequate  when  used  to  represent 
a  small  region  of  a  structure,  but  are  insufficient  when 
behavior  of  the  structure  depends  on  a  small  number  of 
such  elements  (52) .  in  order  to  provide  flexibility  in  an 
analysis  of  concrete  masonry  panels  such  that  it  can  go 
beyond  the  early  cracking  stages,  elements  which  can  be 
adapted  to  general  shapes  should  be  employed.   The  best 
available  element  of  this  kind  is  the  isoparametric 
element  (55) (56) .   Here  the  geometry  and  displacements  of 
the  element  are  described  in  terms  of  the  same  interpola- 
tion functions: 


x  =  {NT}  {Xn}  (7.1) 


u  =  {NT}  {6}  (7.2) 


where  'x'  is  the  x-coordinate  of  any  point  and  'X  '  that 
of  node  'n'  in  an  element,  'u'  is  the  displacement  in 
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x-direction  at  any  point,  and  '6'  represents  the  nodal 
displacements . 

Thus,  the  shape  function  is  identical  to  those  used 
to  generate  the  element.   Equations  of  the  form  of  Eq.  (7.1) 
may  be  considered  as  transformations  between  Cartesian 
coordinates  and  curvilinear  coordinates  if  [N]  is  a  matrix 
of  shape  functions  of  order  two  or  more.   These  equations 
can  be  interpreted  as  the  mapping  of  a  straight-sided  ele- 
ment in  local  coordinates  into  a  curved-sided  element  in 
global  Cartesian  system,  provided  that  the  transformation 
is  unique  in  that  there  is  one-to-one  correspondence  be- 
tween points  in  the  two  coordinate  systems. 

Any  shape  function  chosen,  for  proper  convergence, 
should  provide  for  rigid  body  displacements,  constant 
strain  conditions,  and  finite  strains  at  the  interface 
between  elements.   Zienkiewicz  and  Irons  (56)  have  shown 
how,  for  a  given  number  of  total  degrees  of  freedom,  the 
representation  improves  when  elements  with  progressively 
larger  degrees  of  freedom  are  used  (Figure  7.1).   Accord- 
ingly, a  quadratic  rectangular  element  (8  nodes)  was  chosen 
for  this  investigation  (Figure  7.2).   Satisfying  the  con- 
vergence criteria  mentioned  above,  the  shape  functions 
for  this  element  are: 

Ni  =  1/4  (1  +  £q)  (1  +  m±)     (UL   +   rir)±   -  1) 
for  corner  nodes,  and, 
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Fig.  7.1.  Solution  of  a  cantilever  beam  by 
elements  with  varying  degrees  of 
freedom. 


Fig.  7.2.   Quadratic  rectangular  element, 


142 
KL    =  0:   N.  =  1/2  (1  -  i2)     (1  +  nr^) 

ni  =  0:   N.  =  1/2  (1  +  sq)  (1  -  n2) 

for  midside  nodes. 

7. 2   Formulation  of  Element  Stiffness  Matrix: 

a)   Displacement  function:   in  the  plane  stress  case, 
the  displacements  at  any  point  within  an  element  defined  by 
nodes  1,  2,  3...  can  be  written  as: 


6. 
{p    =    {NL,  N2,  N3...}  (6X}  (7.3) 


where, 


ul 


b)   Strains:   the  strains  at  any  point  can  be 
written  as: 

{e}  =  {B}  {6}e  (7.4) 

For  the  plane  stress  case, 


(7.5) 


9u  |  9u 
9y    9x 
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With  the  shape  functions  [N]  already  obtained,  the 
matrix  [B]  may  easily  be  computed. 

c)   Stresses:   considering  the  initial  strains  e   and 
residual  stresses  a     and  assuming  general  elastic  behavior, 
the  relationship  between  stress  and  strain  will  be  of  the 
form, 


{a}  =  [D]  ({e}  -  ieQ})    +  {aQ} 


(7.6) 


where  [D]  is  the  elasticity  matrix  containing  the 
appropriate  material  properties.   For  the  plane  stress  case, 


{a} 

=    {a*    } 
xy 

1 

V 

0 

[D] 

E 

V 

1 

0 

1   -    v2 

0 

0 

1-v 

(7.7) 


(7.8) 


d)   Equivalent  nodal  forces:   if  a  typical  element 
'e'  is  examined,  the  forces  at  the  nodes  are  uniquely 
defined  by: 

1)   1)   the  displacements  of  those  nodes  (6)  , 

2)  the  distributed  loading  acting  on  the 
element  (p) , 

3)  the  initial  strain  (e  ) ,  and 

4)  the  initial  stress  (a  )  . 

o 
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Assuming  elastic  behavior,  the  characteristic  relationship 
will  always  be  of  the  form, 


{F}e  =  [k]8  {5}e  +  {F}^  +  {F}®  +  {F}^  (7.9) 

oo 


where, 


[k]   represents  the  element  stiffness, 
{F}   the  nodal  forces  required  to  balance  any 
distributed  loads  acting  on  the  element, 

{F}    the  nodal  forces  required  to  balance  any  initial 


strain, 

{F}ff   the  nodal  forces  required  to  balance  any  initial 
o 
stress,  and 

[k]   {6}   forces  induced  through  the  displacement  of  the 
nodes. 


Thus, 

'F 


[F1 


(7.10) 


define  the  equivalent  nodal  forces.   Each  of  the  forces 
(f.)  must  contain  the  same  number  of  components  as  the 
corresponding  nodal  displacement  {6.}  and  be  ordered  in 
the  same  direction. 
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For  the  plane  stress  case, 

in  which  components  U  and  V  correspond  to  the  directions 
of  '  u'  and  'v'  displacements.   X  and  Y  are  the  body  force 
components. 

To  make  the  nodal  forces  statically  equivalent  to 
the  actual  boundary  stresses  and  distributed  loads,  an 
arbitrary  virtual  nodal  displacement  can  be  imposed  and 
the  external  and  internal  work  done  by  the  various  forces 
and  stresses  during  the  displacement  can  be  equated.   Thus, 
we  get, 


{F}e   =    (/    {B}i[D]{B}d(vol)){5}e  -    / {B}i[D]{eo)d(vol) 


+  /  {B}T{aQ}d(vol)  -  /  {N}T{p}d(vol)      (7.11) 


e)   Evaluation  of  element  matrices:   an  evaluation  of 

T 
/  {B}  [D]{B}d(vol)  in  global  coordinates  is  needed  to 

compute  the  element  stiffness.   For  plane  problems, 
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Consider  the  set  of  local  coordinates  and  corresponding 
set  of  global  coordinates  x  and  y, 


3N.    3N.  a     3N.  a 
i  _  1  dx        i  _3y_ 

3£     3x   3£    3y   3? 


9n     3x   3n    3y   3ri 


Thus, 


3N. 

l 

H 

= 

3N. 

l 

3n 

3x  3y 

9?  9£ 

3x  3y_ 

3n  3n 


3N. 


"3~x~ 

3N. 


9y 


9N. 

l 

=  [J] 

9x 

9N. 

l 

9y 

3N. 


3x 
3Ni 


9Y 


=  [J] 


-1 


3N. 

l 


9C 
3N. 


9n 


To  evaluate  [J] 


{x}  =  {NT}{X;.}  ;   {y}  =  {NTHYi} 


3„    3N.        „     3N. 
3x  _  l        _3y_      i 

3£  ~  9^   Xi  ;   9?  "  "3lT  Yi 
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[J]  = 


3N    3N 

2N±       9N2 

9ri    dr\ 


Xl   yl 


(7.12) 


Also,  dx  dy  equals  det  [J]  d£  dn . 

After  finding  the  inverse  of  [J] ,  the  matrix  {B}  could 
be  evaluated.   Numerical  integration  using  the  Gauss 
Quadrature  Formula  is  adopted,  wherein, 

r+1        h 

f  (x)  dx  =  I    H.  f(3j)  (7.13) 

-1  1 

where  H.  are  the  weighting  coefficients,  f(a.)  the  value 
l  3 

of  the  function  at  specified  points  a.,  and  n  the  number 
of  Gauss  points  used. 

7. 3   Material  Properties  and  Failure  Criteria: 

7.3.1  Uniaxial  Stress-Strain  Curves:  the  uniaxial  stress- 
strain  relationships  considered  for  concrete  block,  mortar, 
grout,  and  reinforcing  steel  are  shown  in  Figure  7.3. 

7.3.2  Representation  of  Properties  of  Masonry  Element: 
the  uncracked  masonry  element  is  assumed  to  be  isotropic 
and  the  material  properties  are  represented  by  matrix  [D] : 


[D] 


1 

V 

0 

E 

V 
0 

1 

0 

0 

1-v 
2 

1    -    V2 

(7.14) 
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Fie.  7.3.   Uniaxial  stress-strain  curves. 
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In  the  discretization  of  a  nongrouted  concrete  masonry 
element  for  analysis  element  boundaries  were  made  to  coin- 
cide with  the  edges  of  the  mortar  joints  so  that  the  finite 
element  would  always  represent  a  single  material.   However, 
even  if  the  same  discretization  is  adopted,  every  element 
would  contain  at  least  two  different  materials  in  the  case 
of  a  grouted  masonry  element.   Since  the  problem  could  only 
be  handled  using  a  three-dimensional  analysis,  an  assumption 
has  to  be  made  regarding  the  value  of  E  for  such  elements. 
For  a  plane  problem  this  could  be  approximately  handled  by 
assuming  the  axial  rigidity  of  the  composite  element  as  the 
sum  of  the  axial  ridigities  of  its  components.   No  such 
modification  is  made  for  adjusting  the  Poisson  ratio.   A 
constant  value  is  adopted.   Referring  to  Figure  7.4,  the 
equivalent  modulus,  E   ,  for  the  composite  element  is 
obtained  as, 


Eeq  (at)  =  AG  EG  +  AB  EB  +  Ag  Eg  (7.15) 


where  AQ,  AR,  and  Ag  denote  the  areas  of  grout,  block,  and 
steel  and  E~ ,  ED ,  and  E0  denote  the  respective  moduli.   Per- 
feet  bond  is  assumed  between  steel  and  grout. 

In  order  to  calculate  the  stresses  in  the  components, 
the  following  procedure  is  adopted.     Since  the  strains 
are  the  same  for  each  of  the  components,  the  stresses  in  a 
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Fig.  7.4.   Components  of  a  reinforced 
grouted  element. 


Fie.  7.5.   Rotation  of  coordinates, 


particular  component,  e.g.,  steel,  are  given  by: 
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(7.16) 


In  the  case  of  steel  reinforcement,  the  axial  stresses 
have  to  be  determined  so  as  to  detect  yielding.   Equation 
7.16  can  be  used,  and  the  stress  corresponding  to  the 
direction  of  the  steel  reinforcement  computed.   If,  however, 
the  reinforcement  were  oriented  at  an  angle  to  the  x-y  frame, 
then  using  a  strain  transformation  matrix  [T  ] ,  the  strains 
at  the  new  orientation  could  be  determined.   The  positive 
directions  of  stresses  and  angle  are  indicated  in  Figure  7.5. 

The  strains  corresponding  to  the  new  coordinate  axes 
are  (11)  : 


(      *    1 
e 

X 

P 

•      = 

Y         * 

UXy      J 

cos  cp  sin i  sin<p  cose 
sin2cp  cos24>  -sincp  cost 
■sin2tp      sin2tp      cos2<p 


t,Y 


xy 


(7.17) 


or, 


UM  =  [T  ]{e] 


(7.18) 


where  [T  ]  is  the  strain  transformation  matrix, 
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Now  Equation  (7.16)  is  used  to  compute  the  stress 
in  steel. 

7.3.3   Plastic  Yielding  in  Compression:   the  von  Mises 
yield  criterion  and  the  associated  flow  rule  are  considered 
to  be  applicable  for  biaxial  compression  of  the  masonry 
components.   For  the  plane  stress  case, 


F  =  1/2  [(ox  -  oy)2  +  ay2  +  ax2]  +  3?xy2  -  a^     (7.19) 


where  a  equals  the  uniaxial  stress  at  yield.   The  condition 
that  F  is  less  than,  equal  to,  or  greater  than  zero  estab- 
lishes, respectively,  the  condition  where  the  stress  vector 
is  within  the  yield  surface,  on  the  yield  surface,  or  out- 
side it.   The  material  stiffness  matrix  including  these 
plasticity  conditions  was  derived  in  detail  in  references 
(17)  and  (57).   For  incremental  stress  {Aa}  and  incremental 
strains  {Ae}, 


{Aa}  =  [D   ]{Ae} 


(7.20a) 


where  {Aa}  is  the  change  in  the  stress  vector  which  is 
compatible  with  a  change  in  the  strain  vector,  {Ae},  the 
increment  being  in  the  plastic  region.   The  elasto-plastic 
matrix  [D   ]  is, 


^F 


ID„!     ■.     [Dl     ■      IDHjffyH^riD] 


ep- 


3F    iT 


dF 


h^t^'fTOTl 


-1 


(7.20b) 
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No  plasticity  is  accounted  for  if  the  point  under 
investigation  has  already  failed  under  the  biaxial  com- 
pressive strength  criterion  as  described  in  the  next 
section. 

7.3.4   Biaxial  Strength  of  a  Masonry  Component:   Recently, 
various  researchers  (31  -  35)  have  proposed  failure  enve- 
lopes for  biaxial  strength  of  concrete.   The  failure 
envelope  proposed  by  Liu  et  al .  (33  -  35)  for  biaxial 
compression  is  considered  to  be  applicable  in  this  inves- 
tigation because  of  its  simplicity  and  nondependability 
on  knowledge  of  any  empirical  constants.   Their  experiments 
consisted  of  testing  thin  square  plates  (5"  x  5"  x  1/2") 
under  uniaxial  and  biaxial  compression.   Eight  of  their 
experiments  were  performed  on  idealized  models  of  concrete 
and  eighty- four  used  real  concrete.   Principal  variables 
were  ratio  of  principal  stresses,  water-cement  ratio, 
aggregate-cement  ratio,  and  maximum  size  of  aggregate. 
Special  effort  was  made  to  minimize  the  frictional  confine- 
ment introduced  between  the  bearing  plates  and  concrete. 
The  brush  bearing  platens  designed  by  Buyukozturk  et  al. 

(7)  were  used  in  this  investigation.   Based  on  their  ex- 
perimental results,  the  failure  envelope  shown  in  Figure 

7.6  was  proposed  and  may  be  written  as: 


a  <  0.2  op/aQ    =  1  +   ^  _  ±  (7.21) 
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Fig.  7.6.   Liu's  failure  envelope  for  concrete 
under  biaxial  compression. 
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Fig.  7.7.   Assumed  failure  envelope  for  block, 
mortar,  and  grout. 
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1.0  >  a  >  0.2        ap//ao  =  1*2 

5.0  >  a  >  1.0       o   /o  =    1.2/a                (~A^ 

—   —  p'  o                       cont) 

a  >  5.0  o/an  =  J  (1  + 


p7  o    a  v"    1.2a  -  1' 

where  a   equals  ultimate  strength  of  concrete  plate  in 
uniaxial  compression  in  psi,  a      equals  peak  stress  in 
biaxial  compression  in  psi,  and  a  equals  ratio  of 
principal  stress  in  orthogonal  direction  to  principal 
stress  in  direction  considered. 

The  complete  biaxial  strength  envelope  assumed  for 
the  present  investigation  is  shown  in  Figure  7.7  where 
f   equals  tensile  strength  of  the  component.   While  con- 
sidering grouted  masonry  elements,  the  values  of  oq    and 
a   are  assumed  to  be  as  for  the  weaker  of  the  two  compo- 
nents, viz.,  grout  and  mortar  or  grout  and  concrete  block 
as  the  case  may  be. 

7.3.5  Crushing  of  Masonry:   once  the  compressive  strain 
exceeds  a  limiting  strain,  ecu,  of  the  weaker  component 
in  an  element,  then  it  is  assumed  to  crush  at  that  point, 
and  its  stiffness  is  disregarded,  i.e.,  [D]  is  zero. 

7.3.6  Cracking:   the  principal  stress  planes  can  be 
determined  at  any  point  in  an  element  from  the  stress 
vector  at  that  point.   Since  uncracked  masonry  is  assumed 
to  be  isotropic,  the  constitutive  matrix  [D]  is  applicable 
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in  the  principal  stress  coordinate  system.   The  cracking 
strength  of  masonry  is  determined  from  the  biaxial  failure 
envelope  in  Figure  7.7.   If  the  maximum  principal  tension 
exceeds  this  value,  cracking  will  occur  along  a  plane 
normal  to  the  maximum  principal  stress  direction.   It  is 
assumed  that  the  cracked  concrete  can  no  longer  carry  any 
tensile  forces  perpendicular  to  the  crack.   However,  it  was 
assumed  that  the  irregular  surface  of  crack  could  transfer 
a  certain  amount  of  shear  depending  on  the  value  of  3»  the 
shear  retention  factor.   A  value  of  R  =  0.4  is  assumed  in 
this  investigation.   The  cracked  masonry  is  assumed  to  be 
an  orthotropic  material  with  the  principal  axes  of  or- 
thotropy  (X*-Y*) (Figure  7.8)  normal  and  tangent  to  the 
assumed  crack  direction.   Under  those  conditions  the 
material  stiffness  matrix  for  cracked  component  in  the 
local  principal   coordinate  system  is 


D     --S 

cro 


1  -  v 


0    0      0 
0    1      0 

o   o   e(1--v-)- 


(7.22) 


2 
where  6  equals  shear  retention  factor  (0  <_   B  •  1).   Thus, 


{a*}  =  [DcrQ]  {e*}  (7.23) 
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Fig.  7.8.   Cracking  at  an  integration  point, 
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Fig.  7.9.   Bond  strength  criterion. 
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However,  the  rigidity  matrix  in  Equation  (7.22)  must 
be  referred  to  the  x-y  frame  of  reference.   The  stresses 
and  strains  in  the  X*  -  Y*  frame  are  related  to  the  corre- 
sponding values  in  the  X-Y  frame  through  the  relationships 


{a*}  =  [Tff] {a}  (7.24) 

{£*}  =  [Te]{e]  (7.25) 


where  [T  ]  and  [T  ]  are  the  stress  and  strain  transformation 
matrices,  respectively.   Substituting  Equations  (7.24)  and 
(7.25)  into  Equation  (7.23)  yields 


[T  ]  {a}  -  [D    ]  [T  ]  fe} 
a         cro   e 


or, 


{a}  =  [T  1-1[D    ]  [T  ]{e}  (7.26) 

a      cro    e 


As  the  work  done  by  stresses  on  strains  should  be  the  same 
in  the  global  and  local  systems,  it  can  be  shown  that 
[T  ]    equals  [T£]T.   Therefore, 


{o}  =  [T  ]T[D    ]  [T  ]{e}  (7.27) 

e     cro    e 
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and  the  rigidity  matrix  for  the  cracked  point  in  the 
x-y  frame  is : 


[Dcr]    =     [Te]MDcro][Te] 


(7.28) 


Combining  Equations  (7.28),  (7.22),  and  (7.17) 


1  -  v: 
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sym 


+  6cos2asina 


(sin2a) 


+  6  (cos2a) 

(7.29) 


If  the  masonry  element  had  yielded  plastically  before 
it  cracked,  after  cracking  its  stiffness  is  disregarded. 

7.3.7   Bond  Failure  at  the  Interface  between  Mortar  and 
Block:   debonding  of  mortar  joints  may  occur  in  two  modes, 
a)  compression-shear  and  b)  tension-shear. 

The  failure  criteria  assumed  are  similar  to  those 
adopted  by  Buyukozturk  and  Nilson  (8)  for  modelling  the 
interface  between  mortar  and  coarse  aggregate.   As  described 
in  Chapter  3,  the  compression-shear  bond  strength  can  be 
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modelled  by  the  Coulomb-Mohr  failure  theory  (Figure  7.9) 
Debonding  occurs  when, 


5  =  ?lim  =  c  +  ^an  (7-30) 


where  C  equals  shear  stress  at  the  interface,  a   equals 
normal  stress  at  the  interface,  c  equals  adhesion,  and  y 
equals  coefficient  of  friction  between  mortar  and  block. 

It  is  assumed  that  the  adhesion  is  lost  when  the 
strength  of  the  interface  is  exceeded  and  sliding  is 
resisted  thereafter  only  by  friction: 


%   =   VOn  (7.31) 


Once  £, .   is  exceeded,  to  account  for  the  loss  of  adhesion, 
lim 

zero  shear  modulus  is  assigned  to  the  particular  point. 
In  the  local  coordinate  system,  the  constitutive  matrix 
at  the  point  of  occurrence  of  compression-shear  crack  is 
modified  as  follows: 


[D  .  ]  = 
cbo 
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(7.32) 


Equation  (7.32)  predicts  that  no  shear  can  be  transferred 
by  the  debonded  unit.  Since  the  stresses  are  resisted  by 
friction  in  accordance  with  Equation  (7.31),  a  correction 
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factor  must  be  applied.   Friction  stresses  are  calculated 
depending  upon  the  normal  compression  acting  as  given  by 
Equation  (7.31).   These  stresses  were  replaced  by  stati- 
cally equivalent  forces  at  the  corners  of  the  element, 
details  of  which  are  presented  in  the  next  chapter. 

If  the  joints  were  inclined  at  an  angle  a  to  the 
X-axis,  in  order  to  calculate  the  normal  compressive  stress 
acting  on  the  interface  between  mortar  and  block,  the 
stresses  have  to  be  transformed  to  the  local  system  and 
appropriate  stresses  considered  (Figure  7.10)  depending 
on  the  orientation  of  the  joints.   Also  in  this  case,  the 
matrix  [D  ,  ]  has  to  be  transformed  to  the  global  system 
to  obtain  the  global  constitutive  matrix  [D  .  ]  for  the 
debonded  point.   Using  Equation  (7.28), 
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(7.33) 
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Fig.    7.10.      Normal   stresses   acting  on   interface 
of   inclined  mortar   joints. 
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Presently,  there  is  no  information  on  the  tension-shear 
strength  of  the  interface.   However,  the  modulus  of  rupture 
tests  described  in  Chapter  2  provide  values  for  flexural 
tensile  bond  strength  of  the  mortar  joint.   These  values 
are  assumed  as  applicable  to  pure  tension,  and  a  straight 
line  variation  between  the  pure  tension  and  pure  shear 
failure  points  is  assumed  as  criterion  for  predicting  the 
tension-shear  strength  of  the  interface  (Figure  7.9). 

Due  to  the  nature  of  the  tension-shear  crack,  no  shear 
or  normal  stress  can  be  transmitted  across  the  interface 
after  cracking,  although  compression  or  tension  parallel  to 
the  crack  can  be  carried.   Thus,  the  local  constitutive 
matrix  for  mortar  joint  type  (1)  becomes: 
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xy 


Transformed  to  global  coordinates,  the  constitutive  matrix 

is : 
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For  joint  type  (2) ,  the  local  constitutive  matrix  is 
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(7.36) 


Transformed  to  global  coordinates,  the  constitutive 
matrix  is : 
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7.3.8   Yielding  of  Reinforcement:   using  Equation  (7.16), 
the  stress  in  steel  can  be  determined.   If  the  stress  in 
steel  exceeds  the  yield  limit,  then  E   equals  zero  in 
Equation  (7.15),  and  the  equivalent  modulus  is  modified 
for  subsequent  increments. 


CHAPTER  8 
ANALYSIS  AND  EXAMPLES 

8.1  Method  of  Solution;   The  finite  element  method  used 

in  this  analysis  is  of  the  displacement  type  and  is  applied 

in  an  incremental  form.   Initial  elastic  analysis  continues 

until  either  cracking  or  plasticity  occurs.   The  "initial 

stress"  approach  is  adopted  to  redistribute  the  released 

stresses  due  to  cracking  or  plasticity  to  the  surrounding 

elements.   In  this  analysis,  the  initial  stiffness  is  used 

throughout  each  iterative  cycle  as  well  as  through  all  the 

load  increments.   The  method  is  schematically  illustrated 

in  Figure  8.1.   Load  R^    is  applied  and  a  deflection  U,  is 

computed  using  the  initial  elastic  stiffness  K, .   It  was 

found  that  some  points  have  failed  at  this  level.   This 

changes  the  stiffness  matrix  to  K~.   The  ordinate  A,D, 

represents  the  load  released  due  to  failure  of  some  elements 

(Q1) .   Reloading  the  structure  with  Q,  causes  additional 

displacements  AU,  and  also  additional  cracking  reducing  the 

structure  stiffness  to  K-. .   However,  since  the  same  initial 

stiffness  K,  was  used  again,  the  reloaded  path  (dotted 

lines)  is  D-,B,  instead  of  D,A„.   The  unbalanced  load  on  the 

structure  is  B,D0.   The  process  is  continued  until  no 

further  cracking  occurs  and  the  unbalanced  force  is  within 

a  predetermined  limit. 
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Fig.  8.1.   Schematic  diagram  illustrating  the  nonlinear 
analysis  adopted. 
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8.2   Computational  Procedure:   Stresses  and  strains  are 
computed  at  Gaussian  integration  points.   Three  types  of 
arrangements  were  adopted.   For  very  narrow  elements,  the 
stresses  were  computed  only  at  the  three  Gauss  points 
along  its  longer  centroidal  axis.   However,  when  the  aspect 
ratios  were  smaller,  the  stresses  were  calculated  at  all 
the  nine  points.   For  a  typical  load  increment  'i',  the 
following  procedure  is  adopted: 

1)  The  incremental  displacements  A<5.  and  the  total 
displacements  6.  are  computed.   The  incremental  strains 
Ae.  are  also  computed. 

2)  If  the  Gauss  point  has  not  failed  earlier  due  to 
cracking,  crushing,  or  debonding,  incremental  stresses 
Acr.,  total  stresses  a.,  and  the  principal  stresses  and 
their  direction  are  computed. 

3)  The  minimum  principal  strain  at  the  Gauss  point 
is  checked  to  see  whether  it  exceeds  the  ultimate  compres- 
sive strain  of  the  component.   If  it  does,  then  crushing 
occurs  and  the  entire  stress  vector  at  the  point  (o.  .  + 
Aa . )  is  converted  into  equivalent  nodal  forces  by  using: 


{Q-f}    =    f    {B}T   (ai_1  +  A0;L)  d(vol)  (8.1) 


No  further  computations  are  made  at  this  point  in 
subsequent  increments . 
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4)  If  both  principal  stresses  are  compressive,  it  is 
verified  if  the  peak  stress  exceeds  the  biaxial  compressive 
strength  of  the  component.   In  which  case,  the  point  is 
assumed  to  be  crushed  and  all  the  stresses  are  released 

as  in  the  previous  case.   If  not,  the  next  step  is  executed, 

5)  If  the  point  has  not  failed  before  in  other  modes, 
then  it  is  checked  for  compressive  yielding: 

a)  The  vector  (°i_1  +  Aa i)  is  checked  against 
the  von  Mises  yield  criterion  F.   If  F  <  0,  then  the  point 
is  still  in  the  elastic  state  and  step  (8)  is  executed. 

b)  If  F  >  0  for  the  first  time,  the  value  of  F 

at  the  start  of  the  interval  is  computed  and  the  proportion 
of  the  elastic  stress  and  strain  increment  which  has  taken 
place  in  the  range  above  the  yield  point  is  determined  by 
interpolation  so  that  the  stress  vector  is  now  on  the 
yield  surface.   For  the  above  values  of  elastic  stress,  the 
corresponding  elasto-plastic  stress  increment  is  computed 
using  Equation  (7.20a).   If  Aei_11  is  the  elastic  strain 
increment  added  to  £.  ,  to  obtain  the  stress  vector  on  the 
yield  surface,  the  additional  stress  Ao . '  compatible  with 
the  strain  (Ae.  -  Ae.,1)  is  obtained  from  [D   ]: 

AaiX  =  [Dep] (Aei  "  Aei-l1}  (8'2) 


where  [D   ]  is  the  elasto-plastic  matrix. 


169 


The  stress  and  strain  vectors  are  then  updated  to  be: 


a.P  =   o.       *  +  Aa.1  (8.3) 

l      l-l      i 


:.  ,  +  Ae,  (8.4) 

l-l      l 


where  a.,1  equals  the  stress  vector  compatible  with  F  =  0, 
If  yielding  at  the  Gauss  point  has  occurred  during  an 
earlier  iteration  or  load  increment,  the  compatible  Aa^ 
is  computed  directly  as: 


Aa.1  =  [Dep]  A£i  (8.5) 


The  updated  stress  and  strain  vectors  are  then 
evaluated  as : 


o.p  =  a.  .  +  Aa.1  (8.6) 

i     l-l     i 


£.=£..+  Ae.  (8.7) 

l     l-l      i 


As  the  load  increment  is  finite,  it  is  possible  for 
the  stress  calculated  by  Equations  (8.3)  and  (8.6)  to 
exceed  slightly  the  yield  condition.   This  is  tested,  and 
if  excess  found,  stresses  are  scaled  down  to  remain  on 
the  yield  surface  (o^) . 
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The  initial  stresses  (Aa.  -  (a.  -  cr._.))  are  converted 


into  equivalent  nodal  forces  for  the  element 


{Q2}6  =  /  (B}T  (Aai  -  (ai  -  a^) )  d(vol)         (8.8) 


6)   If  both  principal  stresses  are  tensile,  and  if  the 
maximum  principal  stress  does  not  exceed  the  limiting  ten- 
sile strength  of  the  component,  then  step  (8)  is  executed. 
Otherwise,  the  Gauss  point  is  assumed  to  have  cracked 
normal  to  the  direction  of  the  maximum  principal  stress 
during  this  iteration  and  this  crack  direction  is  locked  in. 
Using  the  elasto-plastic  matrix  for  the  cracked  point,  the 
compatible  stress  vector  is  computed  as: 


{Aai}1  =  [Dcr]  {Ae^  (8.9) 


The  initial  stresses  to  be  supported  by  the  equivalent 
nodal  forces  are: 


{Aai)   =  ({Aai)  -  {Ac^1})  (8.10) 


The  corresponding  nodal  forces  are: 


{Q.}e    =    /{BlMAo.}      d(vol)  (8.11) 

3  1 
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7)  In  the  biaxial  tension-compression  case,  the 
appropriate  failure  criteria  are  applied  to  detect  cracking. 
If  none  is  found,  then  step  (8)  is  executed.   Otherwise, 
the  same  procedure  outlined  in  step  (6)  is  followed. 

8)  It  is  verified  whether  the  point  under  consideration 
is  in  a  mortar  joint.   If  not,  step  (9)  is  executed. 

The  stresses  are  transformed  to  coordinate  axes 
corresponding  to  centroidal  axes  of  the  joints.   The  nature 
of  the  stress  normal  to  the  interface  between  mortar  and 
block  determines  the  mode  of  failure  in  case  the  shear 
stress  exceeds  the  limiting  shear  stress  under  the  corre- 
sponding conditions.   Using  the  appropriate  elasto-plastic 
matrix  (D  .  or  D.,  ) ,  compatible  stresses  are  calculated. 

The  initial  stresses  for  the  case  of  no  previous 
shear-compression  debonding  are: 


{Aai}   =  ({Aai>  -  [D^HAe^)  (8.12) 


The  shear  stresses  resisted  by  friction,  ya  ,  are 
added  to  the  shear  stress  component  of  the  initial  stresses, 
where  a      is  the  normal  compressive  stress  on  the  joint. 

The  equivalent  nodal  forces  are  then  calculated  as: 


{Q5)   =  /{B}i{Aai}   d(vol)  (8.13) 


Similar  criteria  are  applied  for  shear-tension  debonding. 


172 


9)   If  the  point  under  investigation  happens  to  be 
an  element  containing  steel  reinforcement,  the  stresses 
are  checked  for  the  yielding  of  reinforcement.   If  yielding 
occurs,  then  the  equivalent  modulus  of  the  element  is 
modified. 

10)   Steps  (1)  through  (9)  are  carried  out  for  each 
Gauss  point  in  an  element  and  for  each  element  in  the 
structure.   All  the  released  element  nodal  forces  are 
stored  in  their  appropriate  global  locations.   Thus, 

{Q}  =    Z     (Q,e  +  Q9e  +  Q.e...)  (8.14) 

All      x  *  J 

Elements 

If  all  the  elements  of  {Q}  are  less  than  a  predetermined 
value,  then  the  next  load  increment  is  applied  and  steps  (1) 
through  (10)  repeated.   If  not,  the  nodal  forces  {Q}  are 
applied  and  steps  (1)  through  (10)  repeated.   This  process 
is  continued  until  all  elements  of  {Q}  are  below  the  pre- 
determined allowable  nodal  load  which  is  considered  to  be 
negligible  compared  to  the  increment  of  load  applied. 

8.3   Computer  Program:   Figure  8.2  shows  the  general  flow 
chart  of  the  program  adopted.   Figure  8.3  shows  the  flow 
diagram  for  the  computational  procedure  outlined  in  the 
previous  section.   A  modified  Gaussian  Elimination  proce- 
dure was  adopted  for  solving  the  equations  (vide  Appendix  A) . 
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8.4   Examples  of  Nonlinear  Analysis:   The  descriptions  of 
several  test  problems  are  now  presented.   The  details  of 
the  data  deck  required  for  each  problem  are  omitted.   This 
information  is  given  for  one  typical  example  in  Appendix  B 
in  the  form  of  a  sample  computer  output. 

8.4.1   Square  Concrete  Panel  under  Diagonal  Compression: 
this  example  was  chosen  to  check  the  program.   A  square 
concrete  panel  tested  by  Blume  and  Proulx  (6)  was  considered 
for  analysis.   The  panel  measured  48"  x  48"  x  6"  and  had 
no  reinforcement.   The  compressive  strength  of  concrete 
(f  ')  used  was  4350  psi.   The  data  concerning  the  modulus 
of  elasticity  and  tensile  strength  were  not  provided.   The 
above  values  were  assumed,  adopting  the  equations  furnished 
in  the  ACI  Standard  318-71.   Assuming  the  unit  weight  of 
concrete  (w)  as  140  lb/cu  ft, 


E  =  w1,5  33/f-7  =  3.6  x  106  psi 
c 


Splitting  tensile  strength  equals  6.7/f  '  =  440  psi. 

The  maximum  experimentally  observed  shear  stress  for 
the  panel  was  413.5  psi  and  the  corresponding  load  168,416 
lbs.   The  mode  of  failure  for  the  panel  was  a  crack  along 
the  loaded  diagonal  (vertical  diagonal) . 

Figure  8.4  shows  the  finite  element  idealization 
adopted  for  the  panel.   The  stresses  were  determined  at  the 
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Fig.    8. A.      Crack  pattern  of    square   concrete  panel  at    195,000    lbs. 
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Gauss  points  (marked  '  .  '  )  .   The  load  was  applied  in 
increments  of  1,000  lbs,  from  190,000  lbs.   No  cracks 
were  detected  till  the  load  reached  195,000  lbs.   The 
first  cracks  corresponding  to  this  load  are  indicated  in 
Figure  8.4.   All  the  cracks  were  parallel  to  the  loaded 
diagonal.   The  forces  released  due  to  the  cracks  were 
checked  to  see  whether  they  were  within  the  predetermined 
value  of  100  lbs.   If  so,  the  next  increment  was  applied. 
The  crack  propagation  is  illustrated  in  Figures  8.4 
through  8.8.   When  the  load  reached  199,000  lbs,  more 
points  started  cracking,  releasing  large  forces  at  the 
nodal  points.   No  convergence  was  obtained  at  the  end  of 
six  cycles,  indicating  failure.   It  is  evident  from 
Figure  8.8  that  the  analysis  predicted  a  band  of  cracking 
region,  unlike  a  single  line  of  crack  along  the  loaded 
diagonal.   The  ultimate  load  predicted  by  the  analysis 
was  18.2%  higher  than  the  experimental  value  obtained  by  Blume 
and  Proulx   (6) .    If  the  value  of  tensile  strength  assumed 
were  reduced  by  10%,  the  load  corresponding  to  the  initial 
crack  would  be  175,500  lbs.   A  10%  variation  in  split 
cylinder  strength  could  be  anticipated.   However,  the 
analysis  predicted  closely  the  crack  propagation  observed 
in  the  experiment. 

The  solution  time  for  this  problem  was  73.35  seconds 
in  IBM  370  for  the  required  storage  space  of  384  k-bytes. 


178 


/ 

/ 

/  / 
/ 

/ 

/ 

/ 

/ 

/    -     CRACK 

INCLINATION    OF     LINE     INDICATES 

DIRECTION    OF    CRACK 


Fig.    8.5.      Crack  pattern  of    square   concrete   panel   at    196.000    lbs, 
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Fig.    8.6.      Crack  pattern  of   square  concrete  panel   at    197,000   lbs, 
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Fig.    8.7.      Crack  pattern  of   square  concrete   panel   at   198,000   lbs, 
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Fig.    8.8.      Crack  pattern  of   square   concrete   panel   at    199,000   lbs. 
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8.4.2  Hollow  Concrete  Masonry  Panels  under  Diagonal 
Compression:   three  hollow  concrete  masonry  panels  were 
considered.   These  panels  were  made  using  the  model  con- 
crete blocks  and  mortars  MS  or  MW.   The  configurations  of 
the  panels  and  their  experimental  failure  modes  are  shown 
in  Figure  8.9.   All  the  panels  failed  by  debonding  of 
mortar  joints.   For  panels  MS  and  SM,  model  mortar  type 
MS  was  used.   For  panel  MW,  model  mortar  type  MW  was  used. 
The  properties  of  the  mortar  and  blocks  are  given  in 
Chapter  2.   The  experimental  values  of  maximum  loads  for 
panels  SM,  MS,  and  MW  were  300  lbs,  600  lbs,  and  350  lbs, 
respectively. 

In  the  finite  element  analysis,  the  values  of  cohesion 
and  coefficient  of  friction  for  the  model  mortars  were 
assumed  to  be  the  same  as  those  given  by  Equations  3.1  and 
3.2,  respectively,  for  mortars  MW  and  MS.   Similarly, 
equations  3.3  and  3.4  were  assumed  for  obtaining  the  modulus 
of  rupture  values.   Since  the  panels  were  hollow  and  mortar 
joints  were  fully  bedded,  an  equivalent  thickness,  corre- 
sponding to  the  net  area  of  the  block  was  adopted  for 
analysis.   This  thickness  was  calculated  as  the  net  area 
of  the  block  divided  by  the  length  of  the  block.   A  value 
of  1.073"  was  assumed  in  the  analysis. 

Figure  8.10  shows  the  finite  element  discretization 
adopted  for  panel  SM.  Load  was  applied  in  increments  of 
20  lbs,  starting  from  250  lbs.   No  failure  was  detected 
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Fig.  8.9.   Experimentally  observed  failure  modes  in 
hollow  concrete  masonry  panels. 
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Fig.  8.10.   Predicted  failure  pattern  for  panel  SM. 
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till  the  load  was  290  lbs.   At  this  load,  point  1  (Figure 
8.10)  failed  by  compression-shear-debonding.   With  the 
nodal  forces  released  due  to  failure  at  point  1,  the 
second  iteration  was  performed.   This  led  to  the  failure 
of  other  points  shown  in  Figure  8.10.   The  next  iteration 
led  to  the  failure  of  a  large  number  of  points  and  no 
convergence  was  obtained,  marking  failure.   In  order  to 
check  this  failure  load,  the  load  increment  was  reduced  to 
10  lbs.   In  this  case,  also,  the  first  failure  occurred 
at  290  lbs  and  no  convergence  was  obtained. 

Even  though  the  predicted  failure  mode  did  not  agree 
well  with  the  experimental  one,  the  predicted  load  of  290 
lbs  was  close  to  the  average  experimental  value  of  300  lbs. 

The  finite  element  discretizations  for  panels  MS  and 
MW  are  shown  in  Figures  8.11  and  8.12,  respectively.   In 
both  cases,  the  load  was  applied  at  20  lb  increments.   As 
in  the  case  of  panel  SM,  no  failure  was  detected  until  a 
particular  load  was  reached  and  subsequent  iterations  at 
the  same  load  level  did  not  converge,  indicating  failure. 
The  failure  patterns  shown  correspond  to  the  first  two 
iterations . 

It  is  evident  that  the  complete  failure  pattern  could 
not  be  predicted  in  either  case.   However,  the  values  of 
failure  load  predicted  were  quite  close.   For  panel  MS, 
the  predicted  load  of  610  lbs  was  nearly  the  same  as  the 
experimental  value  of  600  lbs.   For  panel  MW,  the 
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Fig.  8.11.   Predicted  failure  pattern  for  panel  MS, 
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Fig.    8.12.      Predicted   failure  pattern   for  panel  MW. 


188 


corresponding  values  are  370  lbs  and  3  50  lbs.   The 
solution  time  for  MS  and  MW  were  115  seconds  and  118 
seconds,  respectively. 

8.4.3   Grouted  Masonry  Piers:   this  example  considers 
the  tests  on  model  piers  described  in  Chapter  4.   The  pier 
chosen  was  PR  3.   The  properties  of  grout,  mortar,  spandrel 
concrete,  etc.  are  summarized  in  Chapter  2.   It  can  be 
noted  that  the  tensile  strengths  of  the  grout  and  block 
were  nearly  equal. 

Figure  8.13  shows  the  finite  element  discretization 
adopted  for  the  pier.   The  numbers  indicate  the  codes  for 
the  material  types.   Table  8.1  summarizes  the  material 
properties  adopted  in  the  analysis.   For  tensile  cracking 
at  a  point,  the  maximum  principal  stress  should  exceed 
the  tensile  strength  of  the  weaker  component.   Accordingly, 
the  value  of  o   was  chosen.   The  equivalent  modulus  was 
calculated  as  described  in  Chapter  7. 

A  preliminary  analysis  predicted  failure  at  about 
30%  of  the  experimentally  obtained  failure  load  because 
of  the  debonding  of  few  vertical  and  horizontal  joints. 
Since  the  grouted  element  consists  of  two  materials  in 
each  element,  the  problem  is  really  three  dimensional  in 
character.   In  the  model  tests,  no  debonding  failure  was 
observed  prior  to  diagonal  cracking.   Hence,  it  was 
decided  to  omit  the  debonding  mode  while  analyzing  grouted 
piers . 
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Fig.    8.13.      Finite   element    idealization   for   pier  PR   3. 
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TABLE  8.1 
MATERIAL  PROPERTIES  USED  IN  ANALYSIS  OF  PR  3 


Material 
Code 

E 

psi 

V 

a 

0 

psi 

at 
psi 

ar 
psi 

0 

1.90 

X 

106 

0.2 

656 

200 

4  2 

2 

2.08 

y 

106 

0.2 

760 

205 

3 

2.22 

X 

106 

0.2 

760 

205 

4 

1.93 

X 

106 

0.2 

656 

200 

5 

1.79 

X 

106 

0.2 

656 

200 

6 

4.78 

X 

106 

0.25 

2300 

450 

7 

5.18 

X 

106 

0.25 

2300 

450 

3 

5.10 

X 

106 

0.25 

2300 

450 

9 

2.10 

X 

106 

0.25 
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In  the  analysis,  the  diagonal  compressive  load  was 
applied  in  increments  of  50  lbs.   No  failure  was  observed 
till  the  load  was  3,800  lbs.   At  3,850  lbs,  the  center 
point  of  the  pier  cracked  first  (marked  '9'  in  Figure  8.14). 
The  crack  direction  was  parallel  to  the  loaded  diagonal. 
At  the  same  load  level,  during  the  second  iteration  due 
to  the  loads  released  by  the  first  cracking,  four  more 
points  (marked  '2')  cracked  in  a  direction  parallel  to  the 
loaded  diagonal.   The  third  iteration,  with  loads  due  to 
the  released  forces,  produced  the  subsequent  crack  pattern 
shown  in  Figure  8.10.   Due  to  the  large  magnitude  of  forces 
released  (30,000  lbs  to  100,000  lbs),  no  convergence  was 
obtained.   Thus,  the  collapse  load  was  taken  as  3,850  lbs. 
The  experimental  value  obtained  was  4,740  lbs.   Considering 
the  variability  of  the  phenomenon,  the  predicted  value  can 
be  considered  satisfactory.   The  predicted  failure  mode 
includes  a  band  of  cracks  parallel  to  the  loaded  diagonal 
instead  of  a  few  well-defined  cracks  obtained  in  the  test. 
However,  the  initial  crack  propagation  is  well-predicted. 
The  subsequent  analytical  divergence  led  to  the  irregulari- 
ties in  the  crack  pattern. 

Figure  8.15  shows  the  analytical  and  experimental  load 
deflection  curves  obtained.   The  analytical  solution  pre- 
sented a  much  stiff er  structure  than  the  actual  one.   A 
similar  conclusion  was  drawn  by  Cervenka  and  Gerstle  in 
their  analysis  of  a  reinforced  concrete  shear  panel  (10) . 
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Fig.    8.14.      Predicted  crack  pattern   for  PR  3, 
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Fig.  8.15.   Analytical  and  experimental  load-deflection  curves 
for  grouted  piers. 
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The  solution  time  for  the  analysis  of  PR  3  was  250 
seconds  for  seven  load  increments  of  50  lbs,  from  3,500  lbs. 

Other  piers  considered  for  analysis  were  PR  7  and 
PR  9.   Table  8.2  shows  the  material  properties  used  in  the 
analysis  of  PR  7.   When  the  debonding  failure  mode  was 
considered,  failure  occurred  at  1,000  lbs.   When  it  was 
not  considered,  the  failure  occurred  at  2,500  lbs,  still 
well  below  the  experimental  value  of  5,600  lbs.   The 
failure  criteria,  based  on  the  weaker  component,  apparently 
were  not  satisfactory  because  of  the  two  different  material 
properties  involved  in  each  element.   In  the  absence  of 
data  concerning  the  actual  distribution  of  stresses  in 
grout  and  mortar,  a  trial  run  was  initiated  adopting  an 
equivalent  tensile  strength  on  the  basis  of  the  areas  of 
of  the  components.   These  values  are  shown  in  parentheses 
in  Table  8.2.   The  analysis  was  performed  by  applying  load 
in  increments  of  50  lbs,  starting  from  2,500  lbs.   Complete 
collapse  occurred  at  2,750  lbs,  about  50%  of  the  experi- 
mental value.   The  points  failed  during  the  first  two  iter- 
ations are  shown  in  Figure  8.16. 

The  material  properties  used  in  the  analysis  of  PR  9 
are  shown  in  Table  8.3.   Since,  in  this  case,  the  grout 
was  weaker  than  both  mortar  and  block,  the  tensile  strength 
of  grout  was  the  governing  criterion  for  cracking.   The 
analysis  was  performed  by  applying  load  at  50  lb  increments 
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TABLE  8.2 
MATERIAL  PROPERTIES  USED  IN  ANALYSIS  OF  PR  7 


.  ,  o  a  .  a 

Material  E  o        t        r 

Code  psi       v       psi      psi      psi 


0 

1.25  x  106 

0.2 

485 

135 

30 

2 

2.16  x  106 

0.2 

2688 

260 

(288) 

3 

2.30  x  106 

0.2 

2688 

260 

(288) 

4 

1.65  x  106 

0.2 

485 

135 

(218) 

5 

1.50  x  106 

0.2 

485 

135 

(218) 
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Fig.    8.16.      Predicted  failure  pattern  for  PR  7. 
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TABLE  8.3 
MATERIAL  PROPERTIES  USED  IN  ANALYSIS  OF  PR  9 


a  o,  o 

o        t        r 

psi      psi      psi 


485  135       33 
450       75 
450       75 
450       75 
450       75 


Material 
Code 

E 
psi 

V 

0 

1 

25  x 

106 

0.1 

2 

1 

85  x 

106 

0.1 

3 

2 

00  x 

106 

0.1 

4 

1 

32  x 

106 

0.1 

5 

1 

18  x 

106 

0.1 
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until  complete  collapse  occurred  at  2,000  lbs.   In  this 
case,  also,  the  predicted  load  was  much  lower  than  the 
experimental  value  of  3,350  lbs.   The  points  failed  in 
the  first  two  iterations  are  shown  in  Figure  8.17. 

The  above  results  indicate  that  the  assumed  plane 
stress  condition  seemed  to  be  valid  only  for  the  case 
where  the  tensile  strengths  of  the  component  materials 
have  comparable  values.   In  all  other  cases,  the  manner  in 
which  the  redistribution  of  the  load  occurs  subsequent  to 
the  failure  of  the  weaker  components  seems  to  be  the  in- 
fluencing factor  in  the  analysis.   It  is  hoped  that  this 
can  be  taken  care  of  in  a  three  dimensional  analysis.   The 
analytical  deflection  curves  indicate  that  the  stiffness 
is  nearly  the  same  for  all  the  three  cases  considered. 
Since  the  values  of  the  various  moduli  are  comparable,  this 
seems  to  be  valid.   However,  the  experimental  results  show 
greater  deviation  for  PR  7  and  PR  9. 

8.4.4   Nongrouted  Masonry  Piers:   two  piers,  PR  6  and  PR  19, 
were  considered  for  analysis.   PR  19  had  an  external  pre- 
compression  of  50  psi  applied  normal  to  bedjoints.   For 
both  the  panels,  an  equivalent  thickness  of  1.073"  was 
adopted  in  the  analysis  for  the  hollow  portion. 

PR  6  consisted  of  mortar  joints  made  of  weak  mortar. 
In  the  analytical  solution  when  the  total  load,  applied  in 
increments  of  50  lbs,  reached  700  lbs,  complete  collapse 
occurred.   The  experimental  failure  load  was  1,100  lbs. 
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Fig.    8.17.      Predicted   failure  pattern  for  PR  9. 
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Fig.  8.18.   Predicted  failure  pattern  for  PR  6. 
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Fig.  8.19.   Predicted  failure  pattern  for  PR  19. 
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The  points  failed  during  the  first  two  iterations  are 
plotted  in  Figure  8.18.   It  is  evident  that  neither  the 
failure  load  nor  the  failure  mode  was  predicted  correctly. 

PR  19  had  mortar  joints  made  of  strong  mortar.   An 
analysis  of  PR  19  revealed  that  complete  collapse  occurred 
when  the  load  was  increased  from  1,350  lbs  to  1,400  lbs. 
The  experimental  value  for  this  case  was  3,200  lbs.   The 
failure  points  corresponding  to  the  first  two  iterations 
are  shown  in  Figure  8.19.   In  this  case,  also,  the 
failure  mode  as  well  as  the  load  was  not  predicted  well. 

It  can  be  concluded  that  the  debonding  criterion 
adopted  was  adequate  to  predict  the  failure  load  correctly 
in  the  case  of  panels  with  fewer  mortar  joints;  they  seem 
to  be  inadequate  in  cases  where  the  joints  are  subjected 
to  bending  and  shear,  as  in  the  case  of  piers.   This  could 
be  due  to  the  confinement  of  the  debonded  regions  provided 
by  the  surrounding  blocks  which  effect  could  not  be  modelled 
with  the  present  experimental  data  available. 


CHAPTER  9 
CONCLUSIONS 


9.1   Conclusions;   An  experimental  and  an  analytical  model 
were  developed  to  analyze  the  problem  of  a  nonreinf orced 
concrete  masonry  shear  wall  with  openings. 

It  was  realized  that  for  a  proper  understanding  of  the 
behavior  of  hollow  masonry  walls,  a  knowledge  of  the  strength 
of  mortar  joints  under  combined  stresses  was  necessary.   An 
experimental  investigation  was  initiated  to  study  the  behav- 
ior of  mortar  joints  under  combined  compression,  shear,  and 
bending.   The  test  specimens  were  standard  3-block  prisms 
with  two  3/8"  mortar  joints.   The  test  variables  were  mortar 
composition  and  strength.   Other  factors,  such  as  suction 
rate,  initial  flow,  block  characteristics  and  curing  proce- 
dures, were  standardized  and  held  constant. 

Test  results  indicate  that  both  shear  strength  and 
flexural  tensile  strength  of  mortar  joints  are  appreciably 
increased  linearly  by  precompression  of  the  joints. 
Combined  flexure  and  shear  under  precompression  can  be  repre- 
sented by  a  circular  interaction  diagram.   The  test  data 
indicate  that  the  current  specifications  are  conservative  in 
limiting  stresses  for  shear  and  flexure- tension. 
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Due  to  limitations  of  space,  capacity  of  testing 
machines,  and  cost,  full-scale  tests  could  not  be  made  to 
study  the  shear  strength  of  concrete  masonry  piers.   One- 
fourth-scale  models  were  used  instead.   The  main  variables 
considered  were  the  grout  and  mortar  strengths  for  the 
grouted  piers;  mortar  strengths  and  magnitude  of  precom- 
pression  were  the  main  variables  for  the  nongrouted  piers. 

The  failure  mode  of  the  grouted  piers  under  diagonal 
compression  loading  was  tensile  cracking  along  the  loaded 
diagonal.   It  was  found  that,  within  the  range  of  mortar 
strengths  adopted,  for  the  same  type  of  grout  used,  a  lower 
strength  of  mortar  did  not  appreciably  alter  the  failure 
load  or  mode.   An  equation  for  obtaining  the  diagonal  ten- 
sile strength  of  the  masonry  was  proposed  by  correlating 
the  average  shear  stress  across  the  cross-section  with  the 
square  root  of  the  corresponding  prism  strength.   It  was 
found  that  the  diagonal  tensile  strength  of  the  masonry 
was  approximately  equal  to  3 . 8  times  the  square  root  of  the 
prism  strength. 

In  the  case  of  hollow  piers,  precompression  increased 
the  shear  capacity  of  the  pier  approximately  linearly.   A 
pier  with  weak  mortar  joints  always  failed  by  debonding  of 
the  mortar  joints,  even  under  high  precompression.   However, 
the  introduction  of  stronger  joints  changed  the  failure 
mode  to  diagonal  tensile  cracking.   Within  the  ranges  of 
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strengths  of  mortar  adopted,  the  mortar  strength  does  not 
appreciably  increase  the  shear  capacity.   It  is  proposed 
that  the  shear  strength  of  a  hollow  pier  may  be  predicted 
by  calculating  the  failure  load  as  the  lesser  of  the  load 
corresponding  to  the  circular  interaction  curve  or  the 
load  corresponding  to  the  maximum  principal  stress  equal 
to  the  tensile  strength  of  the  weaker  component. 

A  nonlinear  finite  element  model  was  chosen  for  the 
analytical  investigation.   An  8-node,  rectangular,  iso- 
parametric element  was  used.   Failure  criteria  were  devel- 
oped to  predict  the  debonding  of  mortar  joints  and  yielding, 
cracking  and  crushing  of  mortar,  grout,  and  block. 

The  finite  element  method  used  was  of  the  displacement 
type  and  was  applied  in  an  incremental  form.   Initial 
elastic  analysis  continued  until  either  debonding,  cracking, 
or  plasticity  occurred.   The  "initial  stress"  approach  was 
adopted  to  redistribute  the  released  stresses  due  to  crack- 
ing to  surrounding  elements.   The  initial  stiffness  was 
used  throughout  each  iterative  cycle  as  well  as  through  all 
the  load  increments.   The  complete  structural  behavior 
could  be  followed  by  means  of  nodal  deflections  and  stresses 
at  integration  points  that  were  printed  after  each  load 
increment.   The  process  was  halted  either  when  the  specified 
load  increments  had  all  been  applied  or  when  the  structural 
deterioration  had  become  sufficiently  extensive  (such  as 
when  very  large  deflections  were  generated) . 
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The  validity  of  the  program  was  checked  by  applying 
it  to  the  study  of  a  square  concrete  panel  subjected  to 
diagonal  compression.   The  program  predicted  the  maximum 
load  within  18%  of  the  experimentally  obtained  value.   The 
application  of  the  analytical  model  was  then  considered 
for  concrete  masonry  panels  and  piers.   The  analytical 
model  predicted  the  collapse  loads  of  nongrouted  square 
concrete  masonry  panels  quite  closely,  but  not  the  failure 
patterns. 

In  the  case  of  grouted  piers,  the  analytical  results 
were  not  satisfactory  when  the  debonding  of  the  joints  was 
considered.   The  predicted  loads  were  only  20%  to  30%  of 
the  experimental  values.   Since  the  grouted  element  consists 
of  two  different  materials,  the  problem  is  one  of  three- 
dimensional  character.   Also,  in  the  model  testing,  no 
debonding  failure  was  observed  prior  to  diagonal  cracking. 
Hence,  it  was  decided  to  omit  the  debonding  mode  of  failure 
while  analyzing  grouted  piers. 

The  collapse  load  was  predicted  within  19%  of  the 
experimental  value  when  the  mortar  and  grout  had  comparable 
tensile  strengths.   When  the  difference  in  strengths  of 
mortar  and  grout  was  considerable,  the  predicted  loads  were 
off  by  about  40%.   It  seems  that  the  manner  in  which  the 
redistribution  of  the  load  occurs  subsequent  to  the  failure 
of  the  weaker  component  is  the  influencing  factor  in  the 
analysis. 
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The  analysis  did  not  give  satisfactory  results  for 
hollow  piers  subjected  to  diagonal  compression.   The  pre- 
dicted loads  were  about  50%  of  the  experimentally  obtained 
values.   Thus,  it  can  be  concluded  that  while  the  debonding 
criteria  adopted  were  adequate  to  predict  the  failure  load 
correctly  in  the  case  of  panels  with  fewer  mortar  joints, 
they  seem  to  be  inadequate  in  cases  where  the  joints  are 
subjected  to  bending  and  shear,  as  in  the  case  of  piers. 

9. 2   Recommendations  for  Further  Study: 

1)  The  circular  interaction  curve  proposed  for 
predicting  the  strength  of  mortar  joints  under  combined 
stresses  was  based  on  tests  on  two  types  of  mortar  joints 
under  a  precompression  of  300  psi.   More  tests  should  be 
conducted  with  different  magnitudes  of  precompression  for 
various  types  of  mortars  to  verify  the  validity  of  the 
equation  over  a  wider  range. 

2)  The  entire  program  on  tests  on  model  piers  should 
be  repeated  for  different  pier  configurations  and  strengths 
of  mortar  and  grout. 

3)  The  analytical  investigations  revealed  a  need  for 
more  satisfactory  criteria  for  predicting  the  debonding  of 
mortar  joints.   Hence,  it  is  recommended  that  more  refined 
criteria  should  be  developed  for  predicting  the  tension- 
shear  debonding  mode. 


APPENDIX  A 
ALGORITHM  FOR  SOLUTION  OF  EQUATIONS 

As  described  in  Chapter  6,  an  algorithm  is  required 
to  solve  equations  of  the  form: 

[K]{6}  =  {R} 

where  [K]  is  the  symmetric  positive  definite  global  stiff- 
ness matrix,  {6}  the  unknown  displacements  and  {R}  the 
right-hand  side  (RHS)  column  vector  of  given  data. 

Since  in  any  problem  there  may  be  some  known 
displacements,  it  is  efficient  if  the  knowns  and  unknowns 
are  separated.   In  this  program,  the  unknowns  are  counted 
column  by  column,  while  the  knowns  are  counted  and  tabu- 
lated separately.   The  stiffness  coefficients  for  knowns 
and  unknowns  in  an  element  are  stacked  up  separately. 

The  resulting  global  stiffness  matrix  will  be 
banded  as  shown  on  the  following  page.   (For  simplicity, 
the  matrix  notations  are  not  adopted  for  the  individual 
matrices. ) 
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Dl 

B1T 

Bl 

D2 
B2T 

CI 
B2 

C1T 

D3 

B3 

D4 

T 
B4 

C3 
B4 

B3T 
C3T 

D5 

B5 
D6 
B6T 

C5 
B6 

B51 
C51 

D7 

C7 


Dn 


xl 

Fl 

x2 

F2 

x3 

F3 

x4 

= 

F4 

x5 

F5 

x6 

F6 

x7 

F7 

xn 

Fn 

where: 

Dl  is  a  square  matrix  of  stiffness  coefficients  which 
relates  the  unknown  displacements  of  first  column  among 
themselves. 

Bl  is  a  matrix  of  stiffness  coefficients  which  relates 
the  unknown  displacements  in  column  1  with  those  of  column  2, 

CI  is  a  matrix  of  stiffness  coefficients  which  relates 
the  unknown  displacements  of  column  1  with  those  of  column  3 
and  so  on. 

xl  is  the  vector  of  unknowns  in  column  1... 

Fl  is  the  corresponding  RHS  vector  for  unknowns  in 
column  1 . . . 

It  is  not  economical  to  store  the  entire  matrix  [K] 
before  solving  for  the  unknowns.   Therefore,  the  following 
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procedure  is  adopted: 

The  matrices  Dl,  Bl...  are  built  up  for  each  column. 
Thus,  when  considering  the  first  column,  the  block, 


Column 

j-l 

J 

j+l 

Dl 

Bl 

CI 

D2 

B2 

D3 

is  first  built  up.   The  transposes  need  not  be  stored. 
However,  D3  is  only  partially  built  up;  the  remainder  of 
D3  is  built  up  when  considering  the  next  column.   Now 
triangularization  can  be  performed  for  this  block. 
Consider  the  first  three  equations: 


Dl  xl  +  Bl  x2  +  CI  x3  =  Fl 


(A.l) 


Bl   xl  +  D2  x2  +  B2  x3  =  F2 


(A. 2) 


C1T  xl  +  B2T  x2  +  D3  x3  +  B3  x4  +  C3  x5  =  F3     (A. 3) 


From  (A.l) , 


xl  =  D1-1(F1  -  Bl  x2  -  CI  x3)  =  Dl  1  Fl' 
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Substituting  for  xl  in  (A. 2), 


(D2  -  BlT  Dl  1  Bl)  x2  +  (B2  -  BlT  Dl  -1    CI)  x3 
(F2  -  B1T  Dl"1  Fl) 


or, 


D2*  x2  +  B2*  x3  =  F2' 


x2  =  (D2*)  1     (F2*  -  B2*  x3) 


Substituting  for  xl  in  (A. 3), 


(B2T  -  C1T  Dl  1    Bl)  x2  +  (D3  -  C1T  Dl"1  Cl)  x3  + 
B3  x4  +  C3  x5  =  (F3  -  ClT  Dl"1  Fl) 


or, 


(B2*)T  x2  +  D3*  x3  +  B3  x4  +  C3  x5  =  F3* 


Substituting  for  x2  in  the  above, 


(D3*  -  B2*T  D2-1  B2)  x3  +  B3  x4  +  C3  x5  =  F3^ 
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or, 


D3**  x3  +  B3  x4  +  C3  x5  =  F3**  (A. 4) 

In  the  preceding  equations,  the  (  )*  and  (  ) **  matrices  are 
defined  implicitly  in  terms  of  the  corresponding  (  )  matrices 
and  are  typically  overstored  in  the  same  core  locations. 

Equation  (A. 4)  can  be  used  to  eliminate  x3  from  the 
next  two  equations,  and  so  on  until  the  last  equation  is  of 
the  form, 

Dn  xn  =  Fn 

from  which  xn  is  obtained.   Back  substitution  yields  the 
other  unknowns . 

The  algorithm  described  above  was  originally  developed 

t 
by  Eisenberg  and  Kim. 


t  Eisenberg,  M.  A.  and  Kim,  B.  S.,  "Numerical  Analysis  of 
Plane  Elastic-Plastic  Boundary  Value  Problems:   Theory 
and  Application  to  Single  Crystal  Beam,"  University  of 
Florida,  Engineering  and  Industrial  Experiment  Station 
Technical  Report  No.  5,  ONR  Contract  No.  N0014-68-A-0173 
-0002,  December,  1970. 


APPENDIX  B 
COMPUTER  PROGRAM  OUTPUT 

The  example  chosen  is  the  square  panel  subjected  to 
diagonal  compression.   The  material  properties  and  loading 
conditions  have  already  been  described  in  Chapter  8. 

Figure  B.l  shows  the  finite  element  idealization  of 
the  panel.   The  computer  program  output  for  the  above 
problem  are  given  in  the  following  pages.   The  input  data 
and  results  of  the  last  load  increment  only  are  given. 
To  facilitate  understanding  the  printout,  the  following 
notes  are  given. 

In  Figure  B.l,  the  secondary  rows  and  columns  are 
shown  in  dotted  lines.   The  topological  codes  for  the 
primary  nodes  are  also  marked  in  the  figure.   The  primary 
nodes  define  an  element.   The  secondary  nodes  are  the 
midside  nodes.   The  integration  points  chosen  are  shown 
in  Figure  8.10.   'DELTA1  denotes  the  shear-retention 
factor.   The  boundary  conditions  are:   1)  at  the  loaded 
point  the  displacement  in  the  x-direction  is  zero; 
2)  at  the  support  the  displacements  in  the  x-  and  y- 
directions  are  zero.   NFL  is  the  matrix  indicating  the 
failure  modes  at  the  integration  points.   'Two  (2)' 
denotes  cracking,  '3',  debonding  (shear-compression),  '4', 
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debonding  (shear-tension),  and  '100',  crushing.   Matrix 
NFL  was  initialized  at  the  start  of  the  analysis.   'TV 
denotes  the  vertical  displacement  of  the  loaded  point. 
Since  the  arrays  TU  and  TV  were  not  initialized,  their 
values  beyond  stage  5  are  not  valid. 
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4    TU 


Fig 


B.l.   Finite  element  idealization  for  panel  SM. 
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